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Abstract. An overview of the author's results is given. Property C stands 
for completeness of the set of products of solutions to homogeneous linear 
Sturm-Liouville equations. The inverse problems discussed include the clas- 
sical ones (inverse scattering on a half-line, on the full line, inverse spectral 
problem), inverse scattering problems with incomplete data, for example, in- 
verse scattering on the full line when the reflection coefficient is known but no 
information about bound states and norming constants is available, but it is a 
priori known that the potential vanishes for x < 0, or inverse scattering on a 
half-line when the phase shift of the s-wave is known for all energies, no bound 
states and norming constants are known, but the potential is a priori known 
to be compactly supported. If the potential is compactly supported, then it 
can be uniquely recovered from the knowledge of the Jost function f(k) only, 
or from /'(0, fc), for all k £ A, where A is an arbitrary subset of (0, co) of 
positive Lebcsgue measure. 

Inverse scattering problem for an inhomogeneous Schrodinger equation is 
studied. 

Inverse scattering problem with fixed-energy phase shifts as the data is 
studied. 

Some inverse problems for parabolic and hyperbolic equations are investi- 
gated. 

A detailed analysis of the invertibility of all the steps in the inversion pro- 
cedures for solving the inverse scattering and spectral problems is presented. 

An analysis of the Newton-Sabatier procedure for inversion of fixed-energy 
phase shifts is given. 

Inverse problems with mixed data are investigated. 

Representation formula for the /-function is given and properties of this 
function arc studied. 

Algorithms for finding the scattering data from the /-function, the I- 
function from the scattering data and the potential from the /-function are 
given. 

A characterization of the Weyl solution and a formula for this solution in 
terms of Green's function are obtained. 



1991 Mathematics Subject Classification. 34B25, 35R30, 81F05, 81F15, 73D25. 
Key words and phrases. Property C, inverse scattering, inverse problems, incomplete data, 
fixed-energy phase shifts, /-function, Weyl function, spectral function. 

1 



Table of Contents 

1. Property C for ODE 

2. Applications of property C 

2.1 Uniqueness of the solution to inverse scattering problem with the data 
/-function. 

2.2 Uniqueness of the solution to inverse scattering problem on the half-line. 

2.3 Compactly supported potentials are uniquely determined by the phase 
shift of s-wave. 

2.4 Recovery of q £ Li i(R) from the reflection coefficient alone. 

2.5 Inverse scattering with various data. 

3. Inverse problems on a finite interval with mixed data. 

4. Property C and inverse problems for some PDE. Recovery compactly 
supported potential from the knowledge of f(k) or f'(0 7 k). 

5. Invertibility of the steps in the inversion procedure, in the inverse 
scattering and spectral problems. 

6. Inverse problem for an inhomogeneous Schrodinger equation. 

7. Inverse scattering problem with fixed-energy data. 

7.1 Three-dimensional inverse scattering problem. Property C. Stability esti- 
mates. 

7.2 Approximate inversion of fixed-energy phase shifts. 

8. A uniqueness theorem for inversion of fixed-energy phase shifts. 

9. Discussion of the Newton-Sabatier procedure for recovery of q(r) from 
the fixed-energy phase shifts. 

10. Reduction of some inverse problems to an overdetermined Cauchy 
problem. An iterative method for solving this problem. 

11. Representation of /-function. 

12. Algorithms for finding q(x) from I(k). 

13. Remarks. 

13.1. Representation of the products of solutions to Schrodinger equation 
(1.1). 

13.2. Characterization of Weyl's solutions. 

13.3. Representation of the Weyl solution via the Green function. 
Bibliography 



PROPERTY C FOR ODE AND APPLICATIONS TO INVERSE PROBLEMS. 3 

1. Property C for ODE 

In this paper a review of the author's results is given and some new results are 
included. The bibliography is not complete. Only the papers and books used in 
the work on this paper are mentioned. The contents of this paper are clear from 
the table of contents. 

The results presented in this paper include: 

1. Property C for ordinary differential equations (ODE), that is, theorems about 
completeness of the sets of products of solutions to homogeneous ODE. 

2. Uniqueness theorems for finding the potential 

a) from the /-function (which equals the Weyl function), 

b) from the classical scattering data for the half-axis problem (a new very 
short proof which does not use the Marchenko method) , 

c) from the phase shift of s-wave in the case when the potential q is com- 
pactly supported and no bound states or norming constants are known, 

d) from the reflection coefficient only (when q = for x < xq), 

e) from mixed data: part of the set of eigenvalues and the knowledge of 
q(x) on a part of a finite interval, 

f) from overdetermined Cauchy data, 

g) from part of the fixed-energy phase shifts, 

h) from various type of data which are typical in PDE problems, 

i) from f(k) or /'(0, k) when q is compactly supported, 

j) from the scattering data for a solution to an inhomogeneous Schrodinger 
equation. 

3. Reconstruction algorithms for finding the potential from overdetermined Cauchy 
data, for finding f(k) and f'(0,k) from the scattering data, for finding the 
scattering data from the /-function and the /-function from the scattering 
data. 

4. Properties of the /-function and a representation formula for it. 

5. Stability estimate for the solution of the inverse scattering problem with fixed- 
energy data. Example of two compactly supported real-valued piecewise- 
constant potentials which produce practically the same phase shifts for all 
values of I. 

6. Discussion of the Newton-Sabatier procedure for inversion of the fixed-energy 
phase shifts. Proof of the fact that this procedure cannot recover generic 
potentials, for example, compactly supported potentials. 

7. Detailed analysis and proof of the invertibility of each of the steps in the 
inversion schemes of Marchenko and Gelfand-Levitan. 

8. Representation of the Weyl solution via the Green function and a characteri- 
zation of this solution by its behavior for large complex values of the spectral 
parameter and x running through a compact set. 

Completeness of the set of products of solutions to ODE has been used for inverse 
problems on a finite interval in the works of Borg and Levitan Jl2| , jl3) . 

Completeness of the set of products of solutions to homogeneous partial dif- 
ferential equations (PDE) was introduced and used extensively under the name 
property C in p6[-[pi|, and JT^j. Property C in these works differs essentially from 
the property C defined and used in this paper: while in p6[-[^l|, and ]l7| property 
C means completeness of the set of products of solutions to homogeneous PDE with 
fixed value of the spectral parameter, in this paper we prove and use completeness 
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of the sets of products of solutions to homogeneous ordinary differential equations 
(ODE) with variable values of the spectral parameter. Note that the dimension of 
the null-space of a homogeneous PDE (without boundary conditions) with a fixed 
value of the spectral parameter is infinite, while the dimension of the null-space 
of a homogeneous ODE (without boundary conditions) with a fixed value of the 
spectral parameter is finite. Therefore one cannot have property C for ODE in the 



sense of [^6|-|3fj, and | J17[ , because the set of products of solutions to homogeneous 
ODE with fixed value of the spectral parameter is finite-dimensional. 

In this paper property C for ordinary differential equations is defined, proved 
and used extensively. Earlier papers are p8| and Q . 

Let 

(1.1) lu := u" + k 2 u - q(x)u = 0, x G R = (-oo, oo). 

Assume 



(1.2) ?elu, Li. 



/OO 
(l + \x\) m \q(x)\dx <oo}. 
-oo 



It is known Q, JTtJ that there is a unique solution (the Jost solution) to (1.1) 
with the asymptotics 

(1.3) f(x,k) = e ikx + o(l), rr^+oo. 

We denote f+(x, k) := f(x, k), f-[x, k) := f(x, -k), tel. The function /(0, k) = 
f(k) is called the Jost function. The function f(k) is analytic in C+ := {k : Imk > 
0} and has at most finitely many zeros in C+ which are located at the points 
ikj,kj > 0, 1 < j < J. The numbers — fcj are the negative eigenvalues of the 

selfadjoint operator defined by the differential expression L q :— — -js + q{x) and 
the boundary condition u(0) = in L 2 (M + ), K + = (0,oo). The function f(k) may 
have zero at k = 0. This zero is simple: if /(0) = then /(0) ^ 0, / := 
Let ip and ip be the solutions to (1.1) defined by the conditions 

(1.4) ^(0,fc) = 0, ^'(0,fc) = l; V(0,fc) = l, ^'(0,fe) = 0, 

where ip' :— It is known |Q, p7| , that tp(x,k) and ip(x,k) are even entire 
functions of k of exponential type < \x\. 

Let g±(x,k) be the unique solution to (1.1) with the asymptotics 

(1.5) g±(x,k) = e ±ikx + o(l), x^-oo 
Definition 1.1. Let p(x) £ Li,i(R + ) and assume 

/•DC 

(1.6) / p(x)f 1 (x,k)f 2 (x,k)dx = 0, Vfc>0, 



Jo 

where fj(x,k) is the Jost solution to (1.1) with q{x) = qj(x), j = 1,2. // (1.6) 

implies p(x) — 0, then we say that the pair {Li, L2}, Lj :~ L qj :— —-j^z +Qj(x) has 
property C+. 

Ifp € Li,i(M_) and 

(1.7) / p(x) gi (x,k)g 2 (x,k)dx = Vfc>0, 

J — OO 

implies p(x) — 0, then we say that the pair {Li, L2} has property C_. 
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In (1.7) gj := c j:j+ {x : k). 

Fix an arbitrary b > 0. Assume that 

(1.8) / p(x)ip 1 (x,k)(p 2 (x,k)dx = Vfc > 
Jo 

implies p{x) — 0. Then we say that the pair {L\, L 2 } has property C v and similarly 
Cip is defined, ipj replacing ipj in (1.8). 

Theorem 1.1. The pair {Li,L 2 }, where Lj :— — j^s +Qj( x ), Qj € L 11 (R + ),j = 
1,2, /ias properties C+,C V and Cy. If qj G ii.i(R_), i/ien {Li,L 2 } has property 
C_. 

Proof. Proof can be found in jl8| . We sketch only the idea of the proof of property 
C+. _ 

Using the known formula 

/>OG 

(1.9) fj{x,k) =e lkx + / A^^yje^di/, j = 1,2, 



where Aj(a;, y) is the transformation kernel corresponding to the potential qj(x), 
j = 1,2, see also formula (2.17) below, and substituting (1.9) into (1.6), one gets 
after a change of order of integration a homogeneous Volterra integral equation for 
p(x). Thus p{x) =0. □ 

The reason for taking b < oo in (1.8) is: when one uses the formula 

. i , , . sin(fcx) f x , . sin(fey) , , „ 

(1.10) cp j (x,k) = — )^ + j o K j (x,y)^^dy, j = l,2, 

for the solution tpj to (1.1) (with q — qj) satisfying first two conditions (1.4), 
then the Volterra-type integral equation for p{x) contains integrals over the infinite 
interval (x, oo). In this case the conclusion p(x) — does not hold, in general. 
If, however, the integrals are over a finite interval (x,b), then one can conclude 
p(x) = 0. 

The same argument holds when one proves property Cy,, but formula (1.10) is 
replaced by 



(1.11) k) = cos(fcx) + / Kj{x, y) cos(ky)dy, j' = l,2, 

Jo 

with a different kernel Kj(x,y). 

2. Applications of property C 

2.1. Uniqueness of the solution inverse scattering problem with the data 
I-function. The I-function I(k) is defined by the formula 

This function is equal to the Weyl function m(fc) which is defined as the function 
for which 

(2.2) W(x,k) := ip(x,k) +m(k)<p(x,k) G L 2 (R + ), Imk > 0, 
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where W[x,k) is the Weyl solution, W(0, k) = 1, W'(0,k) = m(k). Note that 
W(x, k) — ^^k) i as follows from formulas (1.3), (2.1) and from formula (2.3) 
which says I{k) = m(k). 
To prove that 

(2.3) I(k) = m(k), 

one argues as follows. If q 6 Li.i(R + ), then there is exactly one, up to a constant 
factor solution to (1.1) belonging to L 2 (R + ) when Imk > 0. 
Since f(x, k) is such a solution, one concludes that 

(2.4) f(x,k) = c(k)[t(}{x,k)+m{k)ip{x,k)], c(k) ^ 0. 
Therefore, 

_ y/(0,fc) + m(fcV(0,fc) _ 
W ~ V(0,*)+m(*M0,A) W ' 

as claimed. 

In sections 11 and 12 the /-function is studied in more detail. 
The inverse problem (IP1) is: 
Given I(k) for all k > 0, find q(x). 

Theorem 2.1. The IP1 has at most one solution. 



Proof. Theorem 2.1 can be proved in several ways. One way [|19[ is to recover the 
spectral function p(X) from I(k), k — As. This is possible since Iml(k) — 
k > 0, and 



(2.5) dp(X) = { ^/(^)l 2 ' 



A > 0, 



Ei=i c ^(A + rfA, A < 0, kj > 0, 

where — kj are the bound states of the Dirichlet operator L q = —-^2 + q(x) in 
L 2 (M+), f(ikj) = 0, (5(A) is the delta- function, and 

Note that ifcj and the number J in (2.5) can be found as the simple poles of I(k) 
in C+ and the number of these poles, and 

(2.7) Cj = -2ikj Res I(k) = 2k j r 3 , 

where itj := ReSfc=fc J I{k), so r 3 ■ — 

It is well known that dp(X) determines q(x) uniquely jl4[], ftn\ . An algorithm for 
recovery of q(x) from dp is known (Gelfand-Levitan) . In [|l9| a characterization of 
the class of /-functions corresponding to potentials in C™ r (R+ ) , m > is given. 
Here we give a very simple new proof of Theorem 2.1 (cf fig] ): 
Assume that q\ and 52 generate the same I(k), that is, h(k) = h{k) := /(AO- 
Subtract from equation (1.1) for fi(x,k) this equation for f2(x,k) and get: 

(2.8) L 1 w=pf 2 , p(x) := qi{x) - q 2 (x), w := fi(x, k) - f 2 (x, k). 
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Multiply (2.8) by f\ and integrate by parts: 

p(x)f 2 (x, k)h{x,k)dx = {w'h - wf[)\™ = (frf 2 - /(/ 2 )L =0 

(2.9) = hh{h(k) - h{k)) = Vk > 0, 

where we have used (1.3) to conclude that at infinity the boundary term vanishes. 
From (2.9) and property C+ (Theorem 1.1) it follows that p(x) = 0. Theorem 2.1 
is proved. 

□ 

2.2. Uniqueness of the solution to inverse scattering problem on the half 
axis. This is a classical problem fllif , |jl7| . The scattering data are 

(2.10) S={S(k), kj, Sj , 1 < j < J} . 
Here 

(2.11) S(k) := 

is the S-matrix, kj > are the same as in section 2.1, and the norming constants 
Sj are the numbers 

(2.12) Sj := — : ^ > 0. 

fii^f'^ikj) 

Note that (2.6) implies 

(2-13) Cj = Sj lf(0,tkj)} 2 . 

Theorem 2.2. Data (2.10) determine q(x) € ii i i(]R+) uniquely. 

Proof. This result is due to Marchenko Jl4| . We give a new short proof based on 
property C (fl§|l). We prove that data (2.10) determine I(k) uniquely, and then 
Theorem 2.2 follows from Theorem 2.1. To determine I(k) we determine f(k) and 
f'(0,k) from data (2.10). 

First, let us prove that data (2.10) determine uniquely f(k). Suppose there are 
two different functions f(k) and h(k) with the same data (2.10). Then 

< 2 -»> Wi'&Y vkEK - 

The left-hand side in (2.14) is analytic in C + since f(k) and h(k) are, and the 
zeros of h(k) in C + are the same as these of /(fc), namely ikj, and they are simple. 
The right-hand side of (2.14) has similar properties in C_. Thus ^ft is an entire 
function which tends to 1 as |fc| — > oo, so, = 1 and f(k) = h(k). The relation 

(2.15) lim f(k) = 1 

follows from the representation 

(2.16) f(k) = l+f A(0, y)e lk ^dy, A(0, y) e L 1 (M + ). 

Jo 

Various estimates for the kernel A(x, y) in the formula 

poo 

(2.17) f(x,k) = e ikx + A(x,y)e ik ydy 
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are given m 
(2.18) 

(2.19) 

(2.20) 



We mention the following: 
\A(x )V )\<ca '' r + !l 



dA{x,y) 1 / x + y 
dx A q \ 2 

dA(x,y) 1 f x + y 
-^y— + l q 



a(x) :-- 

< ca(x)a 

< ca(x)a 



\q(t)\dt, 

x + y s 
2 

x + y 



where c > here and below stands for various estimation constants. 
From (2.17) and (2.18) formula (2.16) follows. 
Thus, we have proved 

(2.21) S^f(k). 
Let us prove 

(2.22) S^f(0,k). 
We use the Wronskian: 

(2.23) /'(0, k)f(-k) - /'(0, -k)f(k) = 2ik, k e R. 



The function /(fc) and therefore f(—k) — f(k), where the overbar stands for 
complex conjugate, we have already uniquely determined from data (2.10). Assume 
there are two functions f'(0, k) and h'(0, k) corresponding to the same data (2.10). 
Let 

(2.24) w(k) :=/' (0,k) - h'{0,k). 

Subtract (2.23) with h'(0, ±fc) in place of /'(0, ±fc) from equation (2.23) and get 
w(k)f(-k)-w(-k)f(k)=0, 

or 

w(k) w(—k) 



(2.25) 



f(k) f(-k) 



Vfc G 



Claim. j^|y is analytic in C+ and vanishes at infinity and yrz$ * s analytic in 
C- and vanishes at infinity. 



If this claim holds, then j|^y 



0, k € C, and therefore tu(fc) = 0, so /'(0, fc) = 

h'(0,k). 

To complete the proof, let us prove the claim. 
From (2.17) one gets: 



(2.26) 

Taking k 

(2.27) 

Thus 
(2.28) 



f(0, k)=ik- A(0, 0) + / A x (0, y)e ik Vdy 



-oo in (2.16), integrating by parts and using (2.20), one gets: 



1 /" 

--/ A y (0,k)e^dy. 



A(Q,0) = - lim ik[f(k) - 1]. 

k — >oo 
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Since f(k) is uniquely determined by data (2.10), so is the constant A(0,0) (by 
formula (2.28)). 

Therefore (2.24) and (2.28) imply: 

(2.29) lim w(k) = 0. 

|fe|->oo, kec + 

It remains to be checked that (2.10) implies 

(2.30) w(ikj) = 0. 

This follows from formula (2.12): if f(k) and Sj are the same, so are /'(0, ikj), and 
w(ikj) — as the difference of equal numbers 

9i k ■ 

h'(0,ikj) = f'(0,ikj) = — 3 



Theorem 2.2 is proved. □ 

In this section we have proved that the scattering data (2.10) determines the I- 
function (2.1) uniquely The converse is also true: implicitly it follows from the 
fact that both sets of data (2.1) and (2.10) determine uniquely the potential and 
are determined by the potential uniquely. A direct proof is given in section 12 
below. 

2.3. Compactly supported potential is uniquely determined by the phase 
shift of s-wave. Consider the inverse scattering on half-line and assume q{x) = 
for x > a > 0, where a > is an arbitrary fixed number. 

The phase shift of s-wave is denoted by S(k) and is defined by the formula 

(2.31) f(k) = \f(k)\e~ iS ^, 
so the S-matrix can be written as 

(2.32) 5 (fc) = ZL*) =e *«(*). 

f(k) 

If q(x) is real-valued, then 

(2.33) 5(-k) = -6(h), k e M, 
and if q € Li^(R + ), then 

(2.34) 8(oo) = 0. 
Note that S-matrix is unitary: 



(2.35) S(-k) = S(k), \S(k)\ = 1 if keR. 
Define index of S(k): 

(2.36) v := indSVk) := — f dlnStk) = —A m argS{k). 
From (2.32), (2.33) and (2.34) one derives a formula for the index: 

v =-A R 5(k) = ~[S(-0) - S(-oo) + S(+oo) - 5(+0)} 

7T 7T 

(2 ' 37) = _2 ,( +0 ) = J- 2J if ^°)^°' 

TT v ' 1-2J-1 if /(0) = 0. 
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Here we have used the formula: 

(2.38) -S(+0) = #{zcros of f(k) in C+} + ij , 

7T 2 

which is the argument principle. Here So ■= if /(0) ^ and Sq := 1 if /(0) = 0. 
The zero of f(k) at k — is called a resonance at zero energy. 
Let us prove the following result : 

Theorem 2.3. 7/g G i(R + ) decays faster than any exponential: \q(x)\ < ce _c ' :c ' 7 ; 
7 > 1, i/ien i/ie data {S(k) Vfc > 0} determines q(x) uniquely. 

Proof. Our proof is new and short. We prove that, if q is compactly supported 
or decays faster than any exponential, e.g. \q(x)\ < ce~ c ' x ' T , 7 > 1, then S(k) 
determines uniquely kj and Sj, and, by Theorem 2.2, q(x) is uniquely determined. 

We give the proof for compactly supported potentials. The proof for the poten- 
tials decaying faster than any exponentials is exactly the same. The crucial point 
is: under both assumptions the Jost function is an entire function of k. 

If q(x) is compactly supported, q(x) — for x > a, then f(k) is an entire function 
of exponential type < 2a, that is |/(fc)| < ce 2a l fc l (jl|, p. 278]). Therefore S(k) 
is meromorphic in C+ (see (2.32)). Therefore the numbers kj, 1 < j < J, can be 
uniquely determined as the only poles of S(k) in C+. One should check that 

(2.39) /(-**,•) 9^0 if /(**,•)= 0. 

This follows from (2.23): if one takes k = ikj and uses f(ikj) = 0, then (2.23) 
yields 

(2.40) /'(0, ikj) f (-ikj) = -2% < 0. 

Thus f(—ikj) 7^ 0. Therefore S(k) determines uniquely the numbers kj and J. 
To determine Sj, note that 

(2-41) Res S(k)=ft^ = \ Sj , 

k=ik } f(l>kj) 1 

as follows from (2.12) and (2.40). Thus the data (2.10) are uniquely determined 
from S(k) if q is compactly supported, and Theorem 2.2 implies Theorem 2.3. □ 

Corollary 2.1. If q(r) G Li ( x(M+) is compactly supported, then the knowledge of 
f(k) on an arbitrary small open subset of R + (or even on an infinite sequence 
k n > 0, k n ^ k m ifm^n, k n — > k as n — > 00, k > 0) determines q(r) uniquely. 
In section 4 we prove a similar result with the data f'(0,k) in place of f(k). 

2.4. Recovery of q G Li i(IR) from the reflection coefficient alone. Consider 
the scattering problem on the full line: u solves (1.1) and 

(2.42) u~t(k)e lkx , x^+00, 

(2.43) u - e lkx + r(k)e~ lkx , x -» -00. 

The coefficients t(k) and r(k) are called the transmission and reflection coeffi- 
cients (see pj and [0). In general r(k) alone cannot determine q(x) uniquely. 
We assume 

(2.44) q(x) = for x < 0, 

and give a short proof, based on property C, of the following: 
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Theorem 2.4. If q e Li j i(R) and (2-44) holds, then r(fc), Vfc > 7 determines 
q(x) uniquely. 

Proof. We claim that r(k) determines uniquely I{k) if (2.44) holds. Thus, Theorem 
2.4 follows from Theorem 2.1. To check the claim, note that u(x, k) = t(k)f(x, fc), 
so 

and use (2.44), (2.43) to get u = e lkx + r(k)e~ ikx for x < 0, so 

u'(0,k) = ik(l-r(k)) 
1 ' ' u(0,k) l + r(k) ' 

From (2.45) and (2.46) the claim follows. Theorem 2.4 is proved. □ 

2.5. Inverse scattering with various data. Consider scattering on the full line 
(1.1), (2.42)-(2.43), assume q(x) = if x £ [0,1], and take as the scattering data 
the function 

(2.47) u(0,k) := u (k), Vfc > 0. 

Theorem 2.5. Data (2-47) determine q(x) uniquely. 

Proof. If q = for x < 0, then u(x, k) = e %kx + r(k)er lkx for x < 0, u(0, k) = 
1 + r(k), so data (2.47) determine r(fc) and, by Theorem 2.4, q(x) is uniquely 
determined. Theorem 2.5 is proved. Of course, this theorem is a particular case of 
Theorem 2.4. □ 



Remark 2.1. Other data can be considered, for example, u'(0,fc) := v(k). Then 
u'(0,k) = ik[l — r(k)}, and again v(k) determines r(k) and, by Theorem 2.4, q(x) 
is uniquely determined. 

However, if the data are given at the right end of the support of the potential, 
the inverse problem is more difficult. For example, if u(l, k) :— u\(k) is given for 
all k > 0, then u\(k) = t{k)e lk , so t(k) is determined by the data uniquely. 

The problem of determining q{x) from t{k) does not seem to have been studied. 
If q{x) > 0, then the selfadjoint operator L q = —-j^z + q(x) in L 2 (R) docs not have 
bound states (negative eigenvalues). 

In this case the relation |r 2 (fc)| + |t 2 (fc)| = 1 allows one to find 

\r{k)\ = yjl- |t(fc)| 2 , fcGM. 

Define 

(2.48) a(k):=ex P {-^ r^m d s\. 

The function (2.48) has no zeros in C + if L q has no bound states. If q is compactly 
supported then r(fc) and t(k) are meromorphic in C. 
Let us note that 

(2.49) f(x, k) := /+(*, fc) = b(k)g-{x, k) + a{k)g + (x, k) 



(2.50) 



g-{x,k) = c(k)f + (x,k) +d(k)f_(x,k). 
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It is known fjj, |j, that 

(2.51) a(-k)=a(k), b(-k) = ljk), k e M, 

(2.52) c(k) = -b(-k), d(k) = a{k), 

(2.53) a(k)=-^[f(x,k),g.(x,k)], b(k) = ^[f+(x,k),g + (x,k)}, 

where [/, g] :— fg' — f'g is the Wronskian, 

(2.54) |a(£0| 2 = l + |&(£0| 2 , 

The function a(k) is analytic in C+. One can prove jl4| p. 288] 

f°° q(x)dx /]\ 

(2.56) o(fc) = 1 - J -°°^ k fc ^> 

and 

(2.57) b(k)=o(^J, k^oo. 

The function r(k) does not allow, in general, an analytic continuation from the 
real axis into the complex plane. However, if q(x) = for x < 0, then b(k) admits 
an analytic continuation from the real axis into C+ and r(k) is meromorhic in C+. 

If q(x) is compactly supported the functions f±(x,k) and g±{x,k) are entire 
functions of k of exponential type, so that r{k) and t(k) are meromorphic in C. 
From (2.54) one can find \b(k)\ since a(k) is found from formula (2.48) (assuming 
no bound states). 

The conclusion is: recovery of a compactly supported potential from the trans- 
mission coefficient is an open problem. 

3. Inverse problems on a finite interval with mixed data 

Consider equation (1.1) on the interval [0,1]. Take some selfadjoint boundary 
conditions, for example: 

(3.1) -u" + q(x)u - \u = 0, 0<£<1; u(0) = u(l) = 0, A = k 2 . 

Assume q e L l [0, l],g = q. Fix b G (0,1) arbitrary. Suppose q(x) is known on 
the interval [b, 1] and the subset {A m („)} of the eigenvalues of the problem (3.1) is 
known, n — 1,2,... where m(n) is a sequence with the property 

(3.2) Z^ = I(i + e n ), e „^o, a > 0. 

n a 

Theorem 3.1. The data {\ m (n), n = 1> 2, q(x), b < x < 1} determine uniquely 
q(x),0 < x < b, if a > 2b. If a = 2b and Y^=i \ £ \ < °°> then the above data 
determine q(x),0 <x<b, uniquely. 
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Proof. First, assume a > 2b. If there are q\ and qi which produce the same data, 
then as above, one gets 

(3.3) 



G(X) := g(k) := / p(x)ipi(x, k)(p 2 (x, k) dx — (ipiw 1 — ip[w) = Up\w' — <p'\w) 
Jo o 

where w := <p\ — tp2, p '■= qi — (ft, k = ^f~\. Thus 

(3.4) g(k) = at k = ±^ \ m{n) := ±k n . 

The function G(A) is an entire function of A of order | (see formula (1.10) with 
k = \/A) , and is an entire even function of k of exponential type < 2b. One has 

2b\Imk\ 

(3.5) WI^TTf 
The indicator of g is defined by the formula 

(3.6) h(9) := h v (9) := TIm" ln |g(re ' )[ , 

7 — >oo T 

where k = re lS . Since |/mfe| = r\ sin^|, one gets from (3.5) and (3.6) the following 
estimate 

(3.7) h{6) < 2&| sin(9|. 

It is known |l0| formula (4.16)] that for any entire function g(k) ^ of expo- 
nential type one has: 

n(r) 1 f 27r 

(3.8) lim-^< — / h g (6)d6, 

r^oo r Zir J Q 

where n(r) is the number of zeros of g(k) in the disk |fc| < r. From (3.7) one gets 

1 f 2n 2b f 27r 4b 

(3.9) — / h g {9)d0<— \wx0\d8= — 



From (3.2) and the known asymptotics of the Dirichlet eigenvalues: 

(3.10) A n = (nn) 2 + c + o(l), n — > oo, c = const, 
one gets for the number of zeros the estimate 

(3.11) n{r)>2 ^ 1 = 2 — [1 + o(l)], r -> oo. 

^[i+o(^)]<»- 
From (3.8), (3.9) and (3.11) it follows that 

(3.12) cr < 2b. 

Therefore, if a > 26, then c/(fc) = 0. If g(k) = then, by property C v (Theorem 
1.1), p(x) = 0. Theorem 3.1 is proved in the case a > 2b. 
Assume now that a — 2b and 

oo 

(3.13) ^|e„|<oo. 

n=l 

We claim that if an entire function G(A) in (3.3) of order 1 vanishes at the points 

2 2 

(3.14) a„=— 5-(l+e n ), 
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and (3.13) holds, then G(A) = 0. If this is proved, then Theorem 3.1 is proved as 
above. 

Let us prove the claim. Define 



(3.15) 

and recall that 
(3.16) $ (A) 



$(A) 
sin (ex's/A) 



A 



ctVA 

Since G(A„) = 0, the function 



■n i-c 

n—1 v 

n fi 



71=1 



A 



Mr) 



(3.17) 



tu(A) 



g(A) 
d>(A) 



is entire, of order < |. Let us use a Phragmen-Lindelof lemma. 

Lemma 3.1. Jl0| , Theorem 1.22] If an entire function w(\) of order < 1 has the 
property sup- 00<y<ao \w(iy)\ < c, rten u>(A) = c. //, in addition w(iy) — > as 
y — > +oo. £/ien w(A) = 0. 

We use this lemma to prove that w(X) = 0. If this is proved then G(A) = and 
Theorem 3.1 proved. 

The function w(X) is entire of order ~ < 1. 
Let us check that 



(3.18) 



sup \w(iy)\ < oo, 

-oo<y<oo 



and that 

(3.19) \w(iy)\ — > as y — > +oo. 

One has, using (3.5), (3.15), (3.16) and taking into account that cr = 2b: 
(3.20) 

' ' 1 + -S 



){iy)\ = 



G(iy) <f> (iy) 



< 



e 26|JmViy| / e a\Imy/vy\ 



(i + M) 



i 



<- 



1 + I2/I* 



n 



K_ 



< 



1 



\y\ 



n a 



|£nl) < 



Cl 



{ni/j,^^ An} 



Here we have used elementary inequalities: 
(3.21) 



1 + a a , 1 + a 

——, < - if a > d > 0; ; < 1 if < a < d, 

1 + d ~ d ~ 1 + d ~ - - > 



with a := d := M^-, and the assumption (3.13). 



We also used the relation: 

sm{a^/vy) 



a a\Imy/vy\ 



as y 



-oo. 



ay/iy 2a\^/Ty\ 

Estimate (3.20) implies (3.18) and (3.19). An estimate similar to (3.20) has been 
used in the literature (see e.g.j|). 

Theorem 3.1 is proved. □ 
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Remark 3.1. Theorem 3.1 yields several results obtained in Q, and an earlier 
result of Hochstadt-Lieberman which says that the knowledge of all the Dirichlet 
eigenvalues and the knowledge of q(x) on [5,1] determine uniquely q(x) on [0, i] . 

In this case b = i, a = 1. 

One can also obtain the classical result of Borg Q and its generalization due to 
Marchenko @: 

Two spectra (with the same boundary conditions on one of the ends of the 
interval and two different boundary conditions on the other end) determine q(x) 
and the boundary conditions uniquely. 

4. Property C and inverse problems for some PDE 
4.1. Consider the problem 

(4.1) u t = u xx - q(x)u, 0<x <l, <>0, 

(4.2) u(ar,0) = 0, 

(4.3) u(0,t) = 0, u{l,t)=a{t). 

Assume the a(t) ^ is compactly supported, a(t) € L 1 (0,oo), q{x) € i 1 [0, 1], 
problem (4.1) - (4.3) is solvable, and one can measure the data 

(4.4) u'(l, t) := u x {\, t) := b(t) Vt > 0. 

The inverse problem (IP1) is: 
Given {a(t), b(t), Vi > 0} find q(x). 

Theorem 4.1. IP1 has at most one solution. 

Proof. Laplace-transform (4.1) - (4.3) to get 

poo 

(4.5) v" - Xv - q{x)v = < x < 1, v := I e~ xt u(x,t) dt, 

Jo 

poo 

(4.6) w(0,A)=0, v(l, A) = A(X) :— / e~ xt a{t)dt. 



(4.7) t/(l, A) = B(X) := / b(t) e - xt dt. 

Jo 

Assume that there are qi(x) and q2(x) which generate the same data {A(X), B(X), 
VA > 0}. Let p(x) := qi(x) — q2(x),w := vi—V2- Subtract from equation (4.5) with 
v = Vx, q = q\, similar equations with v = V2, q = q2, and get 

(4.8) t\w := w" — Xw — qiw = pv2, 

(4.9) iy(0,A) = 0, w(l, A) = w'(l, A) = 0. 

Multiply (4.8) by tpi(x, A), where 4^1 = 0, yi(0, A) = 0, ^(0, A) = 1, integrate 
over [0, 1] and then by parts on the left-hand side, using (4.9). The result is: 

(4.10) / p(x)v 2 (x,X)(pt(x, X)dx = VA > 0. 



Note that tp\{x, A) is an entire function of A. 
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Since a(t) ^ and is compactly supported, the function A(X) is an entire function 
of A, so it has a discrete set of zeros. Therefore V2(x,X) = c(X)ip2(x, A) where 
c(A) ^ for almost all A e K+. 

Property C v (Theorem 1.1) and (4.10) imply p(x) = 0. Theorem 4.1 is proved. 

□ 

Remark 4.1. One can consider different selfadjoint homogeneous boundary condi- 
tions at x = 0, for example, u'(0, t) — or u'(0, t) — hou(0, t) = 0, ho — const > 0. 

A different method of proof of a result similar to Theorem 4.1 can be found in 
l| and in @. In § some extra assumptions are imposed on q(x) and a(t). 



4.2. Consider the problem: 

(4.11) u t t — u xx - q(x)u, x>0, t>0, 

(4.12) u = u t = at t = 0, 

(4.13) u(0, *)=£(*), 
where <5(t) is the delta-function. Assume that 

(4.14) q(x)=0 if i>1, q = q, geL'[0,l]. 
Suppose the data 

(4.15) u(l,t) := a(t) 

are given for alH > 0. 

The inverse problem (IP2) is: 
Given a(t) Vi > 0, find q(x). 

Theorem 4.2. The IP2 has at most one solution. 

Proof. Fourier-transform (4.11)-(4.13), (4.15) to get 

(4.16) v" + k 2 v - q(x)v = 0, x > 0, 

(4.17) v(0,k) = l, 

(4.18) v(l,k) = A(k) := / a(t)e ikt dt, 

Jq 

where 

/>oo 

(4.19) «(a;,jfe) = / u(x,t)e m dt. 

It follows from (4.19) and (4.16) that 

(4.20) v(x,k) = c(k)f(x,k), 

where f(x,k) is the Jost solution to (4.16). From (4.20) and (4.17) one gets 

(4.21) »<*•'> =ilr 

where f(k) = /(0, k). From (4.21) and (4.18) one obtains 
(4.22, m . «^>. 
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From (4.14) and (4.16) one concludes 

(4.23) f(x,k) =e lkx forx> 1. 
From (4.23) and (4.22) one gets 

(4.24) /( fc ) = 5)- 
Thus /(fc) is known for all k > 0. 

Since q(x) is compactly supported, the data {f{k), Vfc = 0} determine q(x) 
uniquely by Theorem 2.3. Theorem 4.2 is proved. 

□ 

Remark 4.2. One can consider other data at x = 1, for example, the data u x (l, i) 
or u x (l,t) + hu(l,t). The argument remains essentially the same. 



However, the argument needs a modification if (4.13) is replaced by another 

j>,fc) 
/'(o,fe)' 



condition, for example, u x (0,t) — 8(t). In this case v(x,k) = tfc k X , and in place 



of f(k) one obtains /'(0, k) Vfc > from the data (4.18). 

The problem of finding a compactly supported q(x) from the data {/'(0, fc) Vfc > 
0} was not studied, to our knowledge. We state the following: 

Claim. The data /'(0, fc) known on an arbitrary small open subset of (0, oo) or 
even on an infinite sequence of distinct positive numbers k n which has a limit point 
k > 0, determines a compactly supported q(r) € L 1 (K + ) uniquely. 

Our approach to this problem is based on formula (2.23). If /'(0, fc) is known 
for all fc > 0, then /' (0,-fc) = f'(0,k) is known for all fc > 0, and (2.23) can be 
considered as the Riemann problem for finding /(fc) and /(— fc) from (2.23) with the 
coefficients /'(0, fc) and /'(0, — fc) known. If q(x) E Li(M. + ) is compactly supported 
then /'(0, fc) is an entire function of fc. Thus the data determine /'(0, fc), for all 
fc > 0. 

We want to prove that (2.23) defines /(fc) uniquely if /'(0, fc) is known for all 
fc > 0. Assume the contrary. Let /(fc) and h(k) be two solutions to (2.23), and 
w:= f - h, w(k) -> as |fc| -> oo, fc <= C+. Then (2.23) implies 

7^ = ™ — ' 

The function /'(0, fc) has at most finitely many zeros in C+. All these zeros 
arc at the points inj, 1 < j < Ji, where —Kj are the negative eigenvalues of the 

Neumann operator L q := — ^ + q(x) in i 2 (K + ), u'(0) = 0. 

Also /'(0, 0) may vanish. From (2.23) one concludes that w(iKj) = if /'(0, inj) — 
0. Indeed, one has w{inj)f'{Q,—iKj) = w(—iKj)f'(0,iKj). If f'(0,iKj) — then 
/'(0, — inj) 7^ as follows from (2.23). Therefore w(iKj) — as claimed, and the 
function ff^^ is analytic in C+ and vanishes at infinity in C+. Similary, the 
right-hand side of (4.25) is analytic in C- and vanishes at infinity in C_. Thus, 
by analytic continuation, jtjqj^ is an entire function which vanishes at infinity and 
therefore vanishes identically. Therefore w(k) = and /(fc) = h(k). Thus, the data 
{/'(0, fc), Vfc > 0} determines uniquely {/(fc), Vfc > 0}. 

Since q{x) is compactly supported, Theorem 2.3 implies that q{x) is uniquely 
determined by the above data. The claim is proved. □ 
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5. Invertibility of the steps in the inversion provedures in the 
inverse scattering and spectral problems 

5.1. Inverse spectral problem. Consider a selfadjoint operator L q in L 2 (R + ) 
generated by the differential expression L q = —4-? + q{x), q(x) G i ; 1 oc (K + ), q{x) = 
q(x), and a selfadjoint boundary condition at x = 0, for example, w(0) = 0. Other 
selfadjoint conditions can be assumed. For instance: u'(0) — h o u(0) 7 h o = const > 
0. 

We assume that q(x) is such that the equation (1.1) with ImX > 0, A := k 2 , has 
exactly one solution which belongs to L 2 (M. + ) (the limit-point case at infinity). 

In this case there is exactly one spectral function p(X) of the selfadjoint operator 
L q . Denote 

,,,, , u sin(V\x) 

(5.1) ip (x,X):= -j= — . 

Let h{x) € Lg(K+), where L§(M + ) denotes the set of L 2 (R + ) functions vanishing 
outside a compact interval (this interval depends on h(x)). Denote 

/•OO 

(5.2) H(X):= h(x)<p (x,X)dx. 

Jo 

Assume that for every h € Lg(M + ) one has: 

/oo 
H 2 {X)alp{X) = 0^h(x) = 0. 

Denote by V the set of nondecreasing functions p(X), of bounded variation, such 
that if pi,P2 € V, v := p\ — p 2 , and 

(5.4) H:={H(X):heC x (R+)}, 
where H(X) is given by (5.2), then 

(5.5) jy H 2 {X)dv{X)=0 Vh e 7i| => v{X) = 0. 

Theorem 5.1. Spectral functions of the operators L q , in the limit- point at infinity 
case, belong to V . 

Proof. Let b > be arbitrary / <E L 2 (0, b), f = if x > b. Suppose 

/oo 
H 2 (X)dv{X) =0 VheH. 
-oo 

Denote by I+V and I+W the transformation operators corresponding to potentials 
<7i and (/2 which generate spctral functions p\ and p 2l v = p\ — p 2 - Then 

(5.7) = (/ + V)<p! = {I + W)<p 2 , 

where V and W are Volterra-type operators. Condition (5.6) implies: 

(5.8) \\(I + V*)f\\ = \\(I + W*)f\\ V/eL 2 (0,6), 

where V* is the adjoint operator and the norm in (5.8) is L 2 (0, 6)-norm. Note that 

(5.9) Vf:= [ X V(x,y)f(y)dy, 

Jo 
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and 

(5.10) V*f = J V(y,s)f(y)dy. 
From (5.8) it follows that 

(5.11) I + V* =U(I + W*), 

where U is a unitary operator in the Hilbert space H — L 2 (0, b). 

If U is unitary and V, W are Volterra operators then (5.11) implies V — W. 

This is proved in Lemma 5.1 below. If V = W then <pi(x, A) = ipz(x, A), therefore 
qi = qi and pi (A) = piiX). Here we have used the assumption about L q being in 
the limit-point at infinity case: this assumption implies that the spectral function 
is uniquely determined by the potential (in the limit-circle case at infinity there 
are many spectral functions corresponding to the given potential) . Thus if q\ = q%, 
then pi(X) = /02(A). Theorem 5.1 is proved. □ 

Lemma 5.1. Assume that U is unitary and V,W are Volterra operators in H = 
L 2 {0,b). Then (5.11) implies V — W . 

Proof. From (5.11) one gets I + V = (I + W)U* and, using U*U = I, one gets 

(5.12) (I + V)(I + V*) = (I + W)(I + W*). 
Denote 

(5.13) {I + V)' 1 = I +V X , (I + W)- 1 = I + Wi, 
where Vi, W\ are Volterra operators. From (5.12) one gets: 

(5.14) (/ + V*)(I + Wt) = (1+ Vl)(I + W), 
or 

(5.15) V* + + V*W? = Vi + W + V X W. 

Since the left-hand side in (5.15) is a Voleterra operator of the type (5.10) while 
the right-hand side is a Volterra operator of the type (5.9), they can be equal only 
if each equals zero: 

(5.16) V* + W{ + V*W{ = 0, 
and 

(5.17) Vi + W + ViW = 0. 
From (5.17) one gets 

(5.18) Vi(I + W) = -W, 

or [(/ + V)- 1 - + W) = -W. Thus (I + V)-\I + W) = I and V = W as 
claimed. Lemma 5.1 is proved. □ 

The inverse spectral problem consists of finding q(x) given p(X). The uniqueness 
of the solution to this problem was proved by Marchenko |l4j while the reconstruc- 
tion algorithm was given by Gelfand and Levitan [|l) (see also ]l7[|). 

Let us prove first the uniqueness theorem of Marchenko following p3] . In this 
theorem there is no need to assume that L q is in the limit-point at infinity case: 
if it is not, the spectral function determines the potential uniquely also, but the 
potential does not determine the spectral function uniquely. 
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Theorem 5.2. The spectral function determines q(x) uniquely. 
Proof. If qi and q2 have the same spectral function p(X) then 

/CO p OO 

\F 1 (\)fdp(\) = / |F 2 (A)| 2 dp(A)H|.g|| 2 
-oo J —oo 

for any / G £ 2 (0, b), 6 < oo, where 

(5.20) F,(A) := f f(x)<p j (x,\)dx, j = 1,2, 

Jo 

the function <Pj(x, A) solves equation (1.1) with g = and /c 2 = A, satisfies first 
two conditions (1.4), and 

(5.21) g:={I + K*)f, 
where I + K is the transformation operator: 

(5.22) tp2 = {I + K)ip! = if! + [ K(x, y) Vl (y, A) dy. 

Jo 

Note that 

(5.23) F 2 (A)= / f{x){I + K)p 1 dx= / g(x)<pi(x, A) dx. 

Jo Jo 

From (5.19) it follows that 

(5.24) ||/|| = ||(7 + K*)f\\ V/eL 2 (0, &):=#. 

Since Range (I + K*) = H, equation (5.24) implies that I + K is unitary (an 
isometry whose range is the whole space H). Thus 

(5.25) I + K = (I + K*)- 1 = I + T*, 

where T* is a Volterra operator of the type (5.10). 

Therefore K = T* and this implies K = T* = 0. Therefore tp\ = tp2 and q\ = q 2 . 
Theorem 5.2 is proved. □ 

Let dpj(X),j = 1,2, be the spectral functions corresponding to the operators L qj . 
Assume that dpi(X) = cdp 2 (X), where c > is a constant. The above argument 
can be used with a minor change to prove that this assumption implies: c = 1 and 
Q\ =52- Indeed, the above assumption implies unitarity of the operator y/c(I + K). 
Therefore c(I + K) = I + T*. Thus c = 1 and K = T* = 0, as in the proof of 
Theorem 5.2. Here we have used a simple claim: 

If bl + Q = 0, where b = const and Q is a linear compact operator in H, then 
b = and Q = 0. 

To prove this claim, take an arbitrary orthonormal basis {u n } of the Hilbert 
space H. Then — > as n — > oo since Q is compact. Note that ||u„|| = 1, 

so b = 1 1 &ti n 1 1 = ||Qw„|| — > as n — > oo. Therefore b — and consequently Q = 0. 
The claim is proved. 

The Gelfand-Levitan (GL) reconstruction procedure is: 

(5.26) p(X) => L(x, y) => K(x, y) => q{x). 
Here 

/•OO 

(5.27) L{x,y) := \ ipo(x, \)<po{y, X)da(X), da = dp - dp Q , 
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(5.28) dp = 



Compare (5.28) and (2.5) and conclude that po is the spectral function corre- 
sponding to the Dirichlet operator £ q — —-^7 + q{%) in -^ 2 (K+) with q(x) = 0. 
The function K(x,y) defines the transformation operator (cf. (1.10)) 

f x sm(x\/~X) 

(5.29) ip(x, A) = <p (x, A) + J K(x,y)ip (y,X)dy, ip := — -j= — , 

where ip solves (1.1) with k 2 = A and satisifies first two conditions (1.4). 

One can prove (see |llj , jl7) ) , that K and L are related by the Gelfand-Levitan 
equation: 

(5.30) K(x,y)+L(x,y)+ f K(x,t)L(t,y) dt = 0, < y < x. 

Jo 

Let us assume that the data p(X) generate the kernel L(x, y) (by formula (5.27)) 
such that equation (5.30) is a Fredholm-type equation in L 2 (0,x) for K(x,y) for 
any fixed x > 0. 

Then, one can prove that assumption (5.3) implies the unique solvability of the 
equation (5.30) for K(x,y) in the space L 2 (0,x). 
Indeed, the homogeneous equation (5.30) 

(5.31) h(y)+ [ L(t,y)h(t)dt = 0, < y < x, 

Jo 

implies h — if (5.3) holds. To sec this, multiply (5.21) by h and integrate over 
(0, x) (assuming without loss of generaity that h — h, since the kernel L(t,y) is 
real- valued) . The result is: 



0=\\h\\ 2 + \H(X)\ 2 d P (X)- \H(X)\ 2 dp (X), 

J —oc J —00 

or, by ParsevaPs equality, 

/oo 
\H{X)\ 2 dp(X) = 0. 
-00 

From (5.32) and (5.2) it follows that h(y) = 0. Therefore, by Fredholm's alternative, 
equation (5.30) is uniquely solvable. 
If K(x,y) is its solution, then 

(5.33) q(x) = 2^g^. 

If K(x, x) is a C m+1 function then q(x) is a C m -function. 

One has to prove that potential (5.22) generate the spectral function p(X) with 
which we started the inversion procedure (5.26). 
We want to prove more, namely the following: 

Theorem 5.3. Each step in the diagram (5.26) is invertible, so 

(5.34) p ^ L ^ K q. 
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Proof. 1) Step p => L is done by formula (5.27). 

Let us prove L =*> p. Assume there are p\ and p2 corresponding to the same 
L(x,y). Then 



(5.35) 0= / Lp {x 1 X)Lp( ) {y 1 \) dv \/x, y € M. 

J —oc 

Therefore 

/oo 
H 2 {\)dv V/i G C£°( R +)- 
-oo 

By Theorem 5.1 relation (5.36) implies f (A) = 0, so pi = /?2- 

2) Step L =4> if is done by solving equation (5.30) for K(x,y). The unique 
solvability of this equation for K(x,y) has been proved below formula (5.30). 

Let up prove K L. From (5.27) one gets 

, c ,m T! \ L(a; + y) - £(z - y) Z" 00 1 - cos(x\/A) 

(5.37) L(x,y) = - » L ( x ) '■= J ^ da(X). 

Let y = a; in (5.30) and write (5.30) as 

(5.38) L(2x)+ [ K(x,t)[L(x + t) - L(t- x)]dt = -2K(x,x), x>0. 



Note that L{— x) = L(x). Thus (5.38) can be written as: 
(5.39) 

L(2x) + / K(x,s - x)L(s)ds - K(x, x - s)L(s)ds = -2K (x, x), x > 0. 
Jx Jo 

This is a Volterra integral equation for L(s). Since it is uniquely solvable, L(s) 
is uniquely recovered from K(x, y) and the step K =^> L is done. 
3) Step K q is done by equation (5.33). 

The converse step q ^ K is done by solving the Goursat problem: 

(5.40) K xx -q{x)K = K yy 0<y<x, 

1 r 

(5.41) K(x,x) = - / q(t)dt, K(x,0) = 0. 



One can prove that any twice diffcrcntiablc solution to (5.30) solves (5.40)-(5.41) 
with q(x) given by (5.33). The Goursat problem (5.40)-(5.41) is known to have a 
unique solution. Problem (5.40)-(5.41) is equivalent to a Volterra equation (flTij, 

0)- 

Namely if £ — x + y, 77 = x - y, K(x, y) := £?(£, 77), then (5.40)-(5.41) take the 
form 

(5-42) B iv = ±q(^)B(£,r,), Bforj) = \ £ q(t) dt, £(£,£) = 0. 
Therefore 

(5.43) n) = \fq Q) di + i jf g B(f, s) d*. 

This Volterra equation is uniquely solvable for 

Theorem 5.3 is proved. □ 
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5.2. Inverse scattering problem on the half-line. This problem consists of 
finding q(x) given the data (2.10). Theorem 2.2 guarantees the uniqueness of the 
solution of this inverse problem in the class Li.i := £i.i(M+) of the potentials. 

The characterization of the scattering data (2.10) is known (|0, jnj), that is, 
necessary and sufficient conditions on S for S to be the scattering data correspond- 
ing to a q(x) £ L\ \. We state the result without proof. A proof can be found in 
|l7[ . A different but equivalent version of the result is given in Q . 

Theorem 5.4. For the data (2.10) to be the scattering data corresponding to a 
q G L\ i it is necessary and sufficient that the following conditions hold: 

(5.44) i) ind S(k) = -k < 0, n = 2 J or k = 2J+1, 

(5.45) ii) kj > 0, sj > 0, 1 < j < J, 



(5.46) in) S(k) = S{-k) = S- 1 {k), S(oo) = 1, fc e M, 

(5.47) iv) \\F(x)\\ L ^ {m+] + \\F(x)\\ L i {R+) + \\xF'(x)\\ L i {M+) < oo. 

Here k = 2 J + 1 if /(0) = and k = 2 J if /(0) ^ 0, and 

1 r°° J 

(5.48) F(x) :=-z- [1 - S{k)]e lkx dk + ^ s 3 e^ x . 

The following estimates are useful (see @|, p.209, @, JH) , P-569 ): 

\F(2x) +A(x,x)\ < c / \q(x)\dx, \F(2x)\ < c \q{x)\dx, 



\F'{2x)-^p-\<c{i 



where c > is a constant. The Marchcnko inversion procedure for finding q(x) 
from S is described by the following diagram 

(5.49) S => F(x) A(x, y) => q(x). 

The step S => F is done by formula (5.48). 

The step F =>■ A is done by solving the Marchenko equation for A(x, y): 



(5.50) A(x,y)+F(x + y)+ A(x,i)F(t + y) dt = 0, y > x > 0. 

J X 

The step A => q is done by the formula 

, , dA(x.x) 

(5.51) q(x) = -2- 



dx 

It is important to check that the potential q{x) obtained by the scheme (5.4-9) 
generates the same data S with which we started the inversion scheme (5.4-9). 
Assuming q £ Lx,i we prove: 

Theorem 5.5. Each step of the diagram (5.49) is invertible: 

(5.52) S & F & A^ q. 

Proof. 1. The step S =>■ F is done by formula (5.48) as we have already mentioned. 
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The step F S is done by finding kj, Sj and J from the asymptotics of the 
function (5.48) as x — > — oo. As a result, one finds the function 

J 

(5.53) F d (x): hj . 

3 = 1 

If F(x) and Fd{x) are known, then the function 

(5.54) Fs(i) := — / [1 - S(k)]e ikx dk 

is known. Now the function S(k) can be found by the formula 

POO 

(5.55) S(k) = 1 - / F s (x)e~ ikx dx. 



So the step F S is done. 

2. The step F =>■ A is done by solving equation (5.50) for A(x, y). This step is 
discussed in the literature in detail, (see |l4], fl7||). If q S (actually a weaker 
condition Jq 00 a;|g(x)|(ia; < oo is used in the half-line scattering theory), then one 
proves that conditions i) - iv) of Theorem 5.4 are satisified, that the operator 



(5.56) Tf:= / F(y + t)f(t)dt, y>x>0, 

is compact in L l {x, oo) and in L 2 (x, oo) for any fixed x > 0, and the homogeneous 
version of equation (5.50): 

(5.57) / + T/ = 0, y>a;>0, 

has only the trivial solution / = for every x > 0. Thus, by the Fredholm 
alternative, equation (5.50) is uniquely solvable in L 2 (x,oo) and in L 1 (x,oo). The 
step F =>• A is done. 

Consider the step A(x,y) -> F(x). Define 

(5.58) m:=\f<<^ 

The function A{y) determines uniquely f(k) by the formula: 

(5.59) f(k) = l+ I™ A( y )e lk ydy, 



and consequently it determines the numbers ikj as the only zeros of f(k) in C_|_, 
the number J of these zeros, and S(k) = ^gi. To find F(x), one has to find 
Sj. Formula (2.12) allows one to calculate Sj if f(k) and f'(0,ikj) are known. 
To find f'(0,ikj), use formula (2.26) and put k = ikj in (2.26). Since A(x,y) 
is known for y > x > 0, formula (2.26) allows one to calculate f'(0,ikj). Thus 
S(k),kj,Sj,l < j < J, are found and F(x) can be calculated by formula (5.48). 
Step A F is done. 

The above argument proves that the knowledge of two functions A(0, y) and 
Ac (0,y) for all y > determines q(x) £ Li t i(M+) uniquely. 

Note that: 

a) we have used the following scheme A => S => F in order to get the implication 
A^ F, 
and 
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b) since {F(x),x > 0} =>• A and A =>■ {F(x), —oo < x < oo}, we have proved 
also the following non-trivial implication {F(x),x > 0} {F(x), — oo < x < oo}. 

3. Step A =>■ g is done by formula (5.51). The converse step q =>■ A is done by 
solving the Goursat problem: 

(5.60) Acx - q{x)A = A yv , y>x>0, 

1 f°° 

(5.61) A(x,x) = - / q(t)dt, A(x, y) — > as x + y — > oo. 



2 

Problem (5.60)-(5.61) is equivalent to a Volterra integral equation for A(x,y) (see 



17j p.253]). 



2 



(5.62) A{x,y) = \l q(t)dt+l ds I dtq(s - t)A(s - t,s + t). 



One can prove that any twice differentiable solution to (5.50) solves (5.60)-(5.61) 
with q(x) given by (5.51). 

A proof can be found in @, fnj and 0. 

Theorem 5.5 is proved. □ 

Remark 5.1. It follows from Theorem 5.5 that the potential obtained by the scheme 
(5.49) generates the scattering data S with which the inversion procedure (5.49) 
started. 

Similarly, Theorem 5.3 shows that the potential obtained by the scheme (5.26) 
generates the spectral function p(X) with which the inversion procedure (5.26) 
started. 

The last conclusion one can obtain only because of the assumption that q(x) is 
such that the limit-point case at infinity is valid. 

If this is not the case then there are many spectral function corresponding to a 
given q(x), so one cannot claim that the p(X) with which we started is the (unique) 
spectral function which is generated by ?(x), it is just one of many such spectral 
functions. 

Remark 5.2. In [^0| the following new equation is derived: 

/>oo 

(5.63) F(y)+A(y)+ A(t)F(t + y)dt = A(-y), -00 < y < 00, 

Jo 

which generalizes the usual equation (5.50) at x = 0: 

/•OO 

(5.64) F(y) + A(y)+ A(t)F(t + y)dt = 0, y > 0. 



^0 

Since A(—y) = for y > (see (5.58)), equation (5.64) follows from (5.63) for 
y > 0. For y = equation (5.64) follows from (5.63) by taking y — > +0 and using 
(5.58). 

Let us prove that equations (5.63) and (5.64) are equivalent. Note that equation 
(5.64) is uniquely soluble if the data (2.10) correspond to a? € £1,1. Since any 
solution of (5.63) in L 1 (R + ) solves (5.64), any solution to (5.63) equals to the unique 
solution A(y) of (5.64) for y > 0. 

Since we are looking for the solution A(y) of (5.63) such that A = for y < 
(see (5.58)) one needs only to check that (5.63) is satisfied by the unique solution 
of (5.64). 
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Lemma 5.2. Equation (5.63) and (5.64) are equivalent in L 1 (R + ). 

Proof. Clearly, every L X (E + ) solution to (5.63) solves (5.64). Let us prove the 
converse. Let A(y) £ L 1 (R + ) solve (5.64). Define 



(5.65) 



/(£):=!+/ A(y)e ik «dy:=l + A{k) 



We wish to prove that A(y) solves equation (5.63). Take the Fourier transform 
of (5.63) in the sense of distributions. From (5.48) one gets 

/oo J 
F{x)e^ x dx = 1 - S{-i) + 2tt^2 sj 8(£ + ik ) 7 
-oo j=1 

and from (5.63) one obtains: 

(5.67) F(0 + 2(0 + A(-0m = M-$. 

Add 1 to both sides of (5.67) and use (5.65) to get 

(5-68) /(0 + *W("0 = /(-0- 

From (5.66) and (5.68) one gets: 

(5.69) 



f(0 



3=1 



3=1 

Equation (5.69) is equivalent to (5.63) since all the transformations which led from 
(5.63) to (5.69) are invertible. Thus, equations (5.63) and (5.69) hold (or fail to 
hold) simultaneously. Equation (5.69) clearly holds because 

3 ■/ 
(5.70) /(-0 S M + ikj) = E »if( ik i) = 0, 

3=1 3=1 

since ikj are zeros of f(k). 

Lemma 5.2 is proved. □ 

The results and proofs in this section are partly new and partly are based on the 
results in |M| and 0. 

6. Inverse problem for an inhomogeneous Schrodinger equation 
Consider the problem 



3.1) 



u" + k 2 u — q(x)u = —S(x), 



-00 < X < CO, 



d\x\ 



iku — > as Ixl — > oo. 



(6.2) 
Assume 

(6.3) q = q, q = for \x\ > 1, g€L x [-l,l]. 
Suppose the data 

(6.4) {«(-l 1 /c),u(l,fc)}v fc >o 
are given. 

The inverse problem (IP6) is: 
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Given data (6.4), findq(x). 
Let us also assume that 

(A): The operator L q = —-§£1 + q(x) in L 2 (M.) has no negative eigenvalues. 

This is so if, for example, q(x) > 0. 

The results of this section are taken from pl| 

Theorem 6.1. If (6.3) and (A) hold then the data (6.4) determine q{x) uniquely. 

Proof. The solution to (6.1)-(6.2) is 

(fWf/ k s X>Q 

(6 - 5) u = M \i 

where / = f+(x,k), g = g-(x,k), g(k) := g-(0,k), f(k) := f(0,k), [f,g] := 
fg'— f'g = —2ika(k), a(k) is defined in (2.53), / is defined in (1.3) and g is defined 
in (1.5). 

The functions 

are the data (6.4). 

Since q — when x ^ [—1, 1], condition (6.2) implies /(l, k) = e tk , so one knows 

(6.7) hx ty)~A ^^M, V fc>0. 

a[k) a[k) 

From (6.7), (2.49) and (2.50) one derives 

(6.8) 

a{k)h x {k) = -b(-k)f(k) + a(k)f(-k) = -b{-k)h 2 {k)a{k) + h 2 {-k)a(-k)a(k), 
and 

(6.9) a(k)h 2 (k) = b(k)a(k)hi(k) + a{k)h l {-k)a(~k). 
From (6.8) and (6.9) one gets 

(6.10) -b(-k)h 2 {k) + a(-k)h 2 (-k) = hi(k), 
and 

(6.11) b{k)hi{k) + a(-A)fti(-A) = h 2 (k). 
Eliminate b(— k) from (6.10) and (6.11) to get 

(6.12) a(k) = m{k)a(-k) + n(k), Vfc e R, 
where 

Problem (6.12) is a Riemann problem for the pair {a(fc), a(— fc)}, the function 
a(k) is analytic in C+ := {k : k G C, im/c > 0} and <z(— fc) is analytic in C_. The 
functions a(k) and o(— fc) tend to one as k tends to infinity in C + and, respectively, 
in C_, see equation (2.55). 

The function a(k) has finitely many simple zeros at the points ikj, 1 < j < J, 
kj > 0, where — fc| are the negative eigenvalues of the operator £ defined by the 
differential expression tu = —u" + q(x)u in L 2 (R). 

The zeros ikj are the only zeros of a(k) in the upper half-plane k. 
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Define 



1 f 



din a(k). 



— oo 



One has 



(6.14) 



ind a = J, 



where J is the number of negative eigenvalues of the operator £, and, using (6.14) 
and (6.15), one gets 

(6.15) ind m{k) = -2[ind hi(k) + ind h 2 {k)\ = -2[ind g(k) + ind f(k) - 2 J]. 

Since £ has no negative eigenvalues by the assumption (A), it follows that J = 0. 

In this case ind f(k) = ind g(k) = (see Lemma 1 below), so ind m(k) = 0, and 
a(k) is uniquely recovered from the data as the solution of (6.12) which tends to one 
at infinity. If a(k) is found, then b(k) is uniquely determined by equation (6.11) 
and so the reflection coefficient r(k) := is found. The reflection coefficient 
determines a compactly supported q(x) uniquely by Theorem 2.4. 

If q(x) is compactly supported, then the reflection coefficient r(k) := is 
meromorphic. Therefore, its values for all k > determine uniquely r(k) in the 
whole complex /c-plane as a meromorphic function. The poles of this function in 
the upper half-plane are the numbers ikj = 1, 2, . . . , J. They determine uniquely 
the numbers kj, 1 < j < J, which are a part of the standard scattering data {r(k), 
kj, Sj, 1 < j < J}, where Sj are the norming constants. 

Note that if a(ikj) = then b(ikj) ^ 0, otherwise equation (2.49) would imply 
f(x,ikj) = in contradiction to (1.3). 

If r(k) is meromorphic, then the norming constants can be calculated by the for- 
mula sj = —i ^^ = —iReSk^ikj r(k), where the dot denotes differentiation with 
respect to k, and Res denotes the residue. So, for compactly supported potential 
the values of r{k) for all k > determine uniquely the standard scatering data, 
that is, the reflection coefficient, the bound states —kj, and the norming constants 
s j, 1 < j < J- These data determine the potential uniquely. 

Theorem 6.1 is proved. □ □ 

Lemma 6.1. If J = then ind f = ind g = 0. 

Proof. We prove ind f = 0. The proof of the equation ind g = is similar. Since 
ind f(k) equals to the number of zeros of f(k) in C + , we have to prove that f(k) 
does not vanish in C+. If f(z) = 0, z e C+, then z = ik, k > 0, and — k 2 is an 
eigenvalue of the operator £ in L 2 (0, oo) with the boundary condition u(0) = 0. 

From the variational principle one can find the negative eigenvalues of the oper- 
ator £ in L 2 (R + ) with the Dirichlet condition at x — as consequitive minima of 
the quadratic functional. The minimal eigenvalue is: 



where H 1 (K+) is the Sobolev space of iJ 1 (R + )-functions satisfying the condition 
u(0) = 0. 



(6.16) 




O 
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On the other hand, if J = 0, then 

/oo 
[u' 2 + q(x)u 2 ]dx := n u u £ i7 1 (M), ||u|| L 2 (R) = 1. 
-oo 

o 

Since any element u of H 1 (K+) can be considered as an element of iJ 1 (R) if 
one extends u to the whole axis by setting u — for x < 0, it follows from the 
variational definitions (6.16) and (6.17) that Ki < ko- Therefore, if J = 0, then 
Ki > and therefore «o > 0. This means that the operator t on L 2 (]R + ) with the 
Dirichlct condition at x = has no negative eigenvalues. Therefore /(fc) does not 
have zeros in C+, if J = 0. Thus J = implies ind f(k) = 0. 

Lemma 6.1 is proved. □ 

The above argument shows that in general 

(6.18) ind f < J and ind g < J, 
so that (6.15) implies 

(6.19) ind m(k) > 0. 

Therefore the Riemann problem (2.17) is always solvable. It is of interest to 
study the case when assumption (A) does not hold. 

7. Inverse scattering problem with fixed energy data 

7.1. Three-dimensional inverse scattering problem. Property C. The scat- 
tering problem in R 3 consists of finding the scattering solution u := u(x, a, k) from 
the equation 

(7.1) [V 2 + k 2 - q(x)] ip = in R 3 
and the radiation condition at infinity: 

(7.2) V-^o + «, ip --=e lka ' x , afS 2 , 



(7.3) lim 



dv 
o\x\ 



ds = 0. 



I\s\-- 

Here fc > is fixed, S 2 is the unit sphere, a G S 2 is given. One can write 

gifcr / ]\ x 

(7.4) v = A(a', a, fc) ho - as r = \x\ — ► oo, — = a . 

r \r J r 

The coefficient A(a! , a, k) is called the scattering amplitude. In principle, it can 
be measured. We consider its values for a' , a e S 2 and a fixed fc > as the 
scattering data. Below we take fc = 1 without loss of generality. 

Assume that 

(7.5) qeQ a :={q:q = q, q = for \x\ > a, q G L p (B a )} , 

where a > is an arbitrary large fixed number, B a = {x : \x\ < a}, p > |. 

It is known (even for much larger class of the potentials q) that problem (7.1)- 
(7.3) has the unique solution. 

Therefore the map q — » A(a' , a) := A q (a', a) is well defined, 

A(a',a) := A{a , a,fc)U=i- 
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(IP7) The inverse scattering problem with fixed-energy data consists of finding 
q(x) G Q a from the scattering data A{a!,a) Va', a G S 2 . 

Uniqueness of the solution to IP7 for q G Q a (with p = 2) was first announced in 



27 1 and proved in |28j by the method, based on property C for pairs of differential 



operators. The essence of this method is briefly explained below. This method was 



introduced by the author |26| and applied to many inverse problems 27 -]31[], 



m 0L 

In [fL7| and |25| a characterization of the fixed-energy scattering amplitudes is 
given. 

Let {Li,L 2 } be two linear formal differential expressions, 

M, 

(7.6) L jU = ^2 a m j{x)d m u(x), x G M™, n > 1, j = 1,2, 

\m\=0 



where 



g\m\ 

— , |m| = miH \-m Tl 



Let 

(7.7) Nj := Nj (D) := {w : LjW = in D C i?"} 

where 13 is an abitrary fixed bounded domain and the equation in (7.7) is under- 
stood in the sense of distributions. 
Suppose that 



(7.8) / f(x) Wl (x)w 2 (x)dx = 0, wjeNj, feL 2 (D), 



D 



where Wj G Nj run through such subsets of Nj, j = 1, 2, that the products u>iu>2 G 
L 2 (D), and / G L 2 (D) is an arbitrary fixed function. 

Definition 7.1. The pair {Li,L2} has property C if (7.8) implies f(x) = 0, that 
is, the set {wi,w 2 }vw j eN j: Wl w 2 eL 2 (D) is complete in L 2 (D). 

In plj a necessary and sufficient condition is found for a pair {L\, L 2 } with 
constant coefficients, a m j(x) = a m j = const, to have property C (see also jl7J). 

In it is proved that the pair {Li,L 2 } with Lj = —V 2 + q{x), qj G Q a , has 
property C. 

The basic idea of the proof of the uniqueness theorem for inverse scattering 



problem with fixed-energy data, introduced in |27|, presented in detail in 28. , and 
developed in |fl7| , [§9)-f3l]], is simple. Assume that there are two potentials, qi 
and q 2 in Q a which generate the same scattering data, that is, A\ — A 2 , where 
Aj :=A qj (a',a),j = l,2. 
We prove that 0, p.67] 

(7.9) -iir(A 1 -A 2 )= [qi(x) - q 2 (x)]ilj 1 (x,a)tp2(x,-a')dx \/a,a'eS 2 , 



where ipj(x,a) is the scattering solution corresponding to qj, j = 1, 2. 
If A\ = A 2 , then (7.9) yields an orthogonality relation: 

(7.10) / p{x)4>\{x, a)ip 2 {x, P) dx — Vq;,/3g5' 2 , p(x) := q\ — q 2 . 

Jd 
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Next we prove |l7], p.45] that 

(7.11) spaa a&s3 {ip(x,a)} is dense in L 2 (D) in Nj(D) H H 2 {D), 
where H m {D) is the Sobolev space. Thus (7.10) implies 

(7.12) / p(x)w 1 {x)w 2 (x)dx = Vuij G N j (D)nH 2 (D). 
Jd 

Finally, by property C for a pair {L\, L 2 }, Lj — — V 2 +<&, % G Q a , one concludes 
from (7.11) that p(x) = 0, i.e. q\ = q 2 . We have obtained 

Theorem 7.1 (Ramm [[||, 0). TTie data A(a',a) Va', a G 5 2 determine q G 
Q a uniquely. 

This is the uniqueness theorem for the solution to (IP7). In fact this theorem is 
proved in [T^| in a stonger form: the data Ai(a',a) and A2(a',a) are asumed to 
be equal not for all a! , a G S 2 but only on a set S 2 x S 2 , where S 2 is an arbitrary 
small open subset of S 2 . 

In |24[] the stability estimates for the solution to (IP7) with noisy data are ob- 
tained and an algorithm for finding such a solution is proposed. 

The noisy data is an arbitrary function A s (a!,ot), not necessarily a scattering 
amplitude, such that 

(7.13) sup \A(a',a)-A e (a',a)\<£. 

a',aeS 2 

Given A s (a' , a), an algorithm for computing a quantity q e is proposed in [ p4[ , such 
that 

(7.14) su P \q e - m \<c^l. 

where c > is a constant depending on the potential but not on e, 

(7.15) q(0 := f e* x q(x)dx. 

JB a 

The constant c in (7.14) can be chosen uniformly for all potentials q G Q a which 
belong to a compact set in L 2 (B a ). 

The right-hand side of (7.14) tends to zero as e — > 0, but very slowly. 

The author thinks that the rate (7.14) cannot be improved for the class Q a , but 
this is not proved. 

However, in jjj an example of two spherically symmetric piecewise-constant po- 
tentials q(r) is constructed such that \q\ — q%\ is of order 1, maximal value of each 
of the potentials q is of order 1, the two potentials are quite different but they 
generate the set of the fixed-energy (k = 1) phase shifts > 3 = 1, 2, 

such that 



(7.16) S^=sf\ 0<£<4, 



< 10" 5 , V£ > 5. 



In this example e ~ 10^ 5 , (In | lne|) 2 ~ 2.59, | lne| ~ 5, so the right-hand side of 
(7.14) is of order 1 if one assumes c to be of order 1. 

Our point is: there are examples in which the left and the right sides of estimate 
(7.14) are of the same order of magnitude. Therefore estimate (7.14) is sharp. 
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7.2. Approximate inversion of fixed-energy phase shifts. Let us recall that 
q(x) = q(r), r — \x\ if and only if 

(7.17) A(a',a,k) = A(a'-a,k). 

It was well known for a long time that if q — q(r) then (7.17) holds. The converse 
was proved relatively recently in |32| (see also [p7| ) . 
If q = q(r) then 

oo 

(7.18) A(a', a)=J2 AgYg{a')Yja), 

1=0 

where Yg are orthonormal in L 2 (S 2 ) spherical harmonics, Yg — Yg m , —£< m < £, 
summation with respect to m is understood in (7.18) but not written for brevity, 
the numbers Ag are related to the phase shifts Sg by the formula 

(7.19) A e = Aire iSt sin(^) (k = 1), 
and the S- matrix is related to A{a' , a, k) by the formula 

S = I-—A. 

If q = q(r), r — \x\, then the scattering solution ip(x, a, k) can be written as 



(7.20) ^ {x ,a,k)=T^i*^^Ye(x°)Wa), x°:= X -. 



1=0 



The function ipe( r > k) solves (uniquely) the equation 

/>oo 

(7.21) ipe(r,k) =u e (kr) - g e {r, p)q{p)^g{p 1 k) dp, 

Jo 

where 



(7.22) u t {kr) := ^ — J e+1/2 (kr), v e := y — N e+1/2 (kr). 

Here Jg and Ng are the Bessel and Neumann functions, and 

{ <Poe(kp)foe(kr) f > 
Foi(k) ' ' - P> 

(7-24) F oe (k) = 

ipse, and fog solve the equation 

(7.25) + feVo^ - ^J 1 ^ = 0, 
and are defined by the conditions 

(7.26) f oe ~ e tkr as r -> +oo, 
so 

(7.27) / ol (fcr)=«e^(Mf + ^), 
where it£ and are defined in (7.22), and 

(7-27') ^(fcr) := 
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In ]3S| ] an approximate method was proposed recently for finding q(r) given 
{h}t=o,i,2,---- I n © numerical results based on this method are described. 

In physics one often assumes q{r) known for r > a and then the data {Se}i=o,i,2,... 
allow one to calculate the data ipi(a),£ — 0, 1, 2, ... , by solving the equation 

£(£ + 11 

(7.28) $ + fa - ^ V / " «(r)^ = 0, r > a 
together with the condition 

(7.29) fa ~ e l5f sin(r- — ' + r -> +oo, fc = 1, 

and assuming q(r) known for r > a. 

Problem (7.28) and (7.29) is a Cauchy problem with Cauchy data at infinity. 
Asymptotic formula (7.29) can be differentiated. If the data ipe(a) are calculated 
and ipi{ r ) fo r r > a is found, then one uses the equation 

(7.30) 1 p i ( r )=4°\r)- f 9i (r, p)q(p)Mp) d P> < r < a, 

Jo 

where gt is given in (7.23), 

/■oo 

(7.31) ipf \r):=u e (r)- g e (r, p)q(p)ipt(p) dp, 

J a 

and ui(r) is given in (7.22). 
Put r = a in (7.30) and get 

(7.32) f g l (a 1 p)Mp)<l(p)dp^4 0) ( a )-'Ma) ~h, ^ = 0,1,2. 
Jo 

The numbers bi are known. If q{p) is small or I is large then following approxi- 
mation is justified: 

(7.33) Mp)~4 0) (p)- 
Therefore, an approximation to equation (7.32) is: 



(7.34) / f l {p)q{p)dp = b l , £ = 0,1,2,..., 

Jo 

where 

(7-35) Mp):=g e (a, P )4°\p). 

The system of functions {fe(p)} is linearly independent. 

Equations (7.34) constitute a moment problem which can be solved numerically 
for q(p) (@ p.209], |0|). 

8. A UNIQUENESS THEOREM FOR INVERSION OF FIXED-ENERGY PHASE SHIFTS. 

From Theorem 7.1 it follows that if q = q(r) £ Q a then the data 

(8.1) {S e } W = 0,1,2,... k=l 

determine q(r) uniquely p8]| . 

Suppose a part of the phase shifts is known (this is the case in practice) . 
What part of the phase shifts is sufficient for the unique recovery of q(r) ? 
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In this section we answer this question following |23[ ] . Define 

(8.2) £:=h:J2j = ™ 

\ etc 

to be any subset of nonnegative integers such that condition (8.2) is satisfied. 
For instance, C = {2£}^ = o,i,2,... or £ = {21+ l}^=o,i,2,... will be admissible. 
Our main result is 

Theorem 8.1 ((2^]). If q(r) €E Q a and L satisfies (8.2) then the set of fixed- energy 
phase shifts {^}v£e£ determine q(r) uniquely. 

Let us outline basic steps of the proof. 

Step 1. Derivation of the orthogonality relation: 

If qi, 52 € Qa generate the same data then p(r) := q\ — q2 satisfies the relation 



(8.3) / •p{r)ipf ) {r)i)f ) {r)dr = Q W e £. 



Here ipP (r) are defined in (7.20) and correspond to q — qj, j — 1,2. Note that 
Ramm's Theorem 7.1 yields the following conclusion: if 

(8.4) / p{r)il)f\r)i)f ) (r)dr = Q W e £, 
J o 

then p(r) — 0, and qi = q%. 

Step 2. Since = c^ip^(r), where Sp are some constants, relation (8.3) is 
equivalent to 

(8.5) / p{r)if [ p (r)ipf ] (r) dr = W G £. 
Jo 

Here w'^(r) is the solution to (7.28) with q = qj, which satisfies the conditions: 

(8-6) ^ = (27Tl)!! +0(r ' +1) ' r ^°' 

and 

(8.7) tpf = sin(r - y + + o(l), r -> +oo, 

where \Fp \ ^ is a certain constant, and <rf^ is the fixed-energy (fc = 1) phase 

shift, which does not depend on j by our assumption: = W € £. 

We want to derive from (8.5) the relation (7.10). Since q(x) — q(r), r = \x\ in 
this section, relation (7.10) is equivalent to relation (8.12) (see below). The rest of 
this section contains this derivation of (8.12). 

We prove existence of the transformation kernel K(r, p), independent of £, which 
sends functions u^{r), defined in (7.22), into (fe(r): 

(8.8) Mr)=u e (r)+ f K(r,p)u e (p)% K(r,0)=0. 

Jo P 

Let 

(8.9) 7 , := y| r QW^ + i), 
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where T(z) is the gamma function, 

(8.10) H{1) := 7 | f p(r)ul(r)dr, 

Jo 

and 

(8.11) h(£):=j* f p{r)^\r) V f\r)dr. 

Jo 

We prove that: 

(8.12) {h(£) = W e £} => {h(£) = W = 0, 1, 2, . . . } . 

If h(£) = 0, W = 0, 1, 2, . . . , then (8.4) holds, and, by theorem 7.1, p(r) = 0. 
Thus, Theorem 8.1 follows. 
The main claims to prove are: 

1) Existence of the representation (8.8) and the estimate 

(8.13) [ \K(r, p) I— < c(r) < oo Vr > 0. 
Jo P 

2) Implication(8.12). 

Representation (8.8) was used in the physical literature (||, |l6|]) but, to our 
knowledge, without any proof. Existence of transformation operators with kernels 
depending on £ was proved in the literature ji3) . For our purposes it is important 
to have K(r, p) independent of £. 

Implication (8.12) will be estblished if one checks that h(£) is a holomorphic 
function of £ in the half-plane J\ + :—{£:£ E C, Re£ > 0} which belongs to 
N-Class (Nevanlinna class). 

Definition 8.1. A function h(£) holomorphic in Y[ + belongs to N-class iff 

' 1 - re lv ~ 



(8.14) sup / ln H 

0<r<l , 



1 + re % f 



dip < oo. 



Here 

, I u if u > 0, 
:= { 

[0 if u < 0. 

The basic result we need in order to prove (8.12) is the following uniqueness theo- 
rem: 

Proposition 8.1. If h(£) belongs to N-class then (8.12) holds. 

Proof. This is an immediate consequence of the following: 

Theorem ( |4l| , p. 334]/- If h(z) is holomorphic in D\ := {z : \z\ < 1, z E C}, 
h(z) is of N-class in D±, that is: 

ln+ \h(re lv )\dp < oo, 

and 

(8.16) h(z n ) = 0, h = 1,2,3,..., 

where 

oo 

(8.17) £(1 - |*„|) = oo, 

h=l 
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then h{z) = 0. 

The function z = y=| maps conformally J\ + onto D\, £ = and if h{l) = 
W G £, then /(z) := /i(fef) is holomorphic in D\, f{zi) — for £ S £ and 
^ : = i^f , and 

1 



(8.18) 



1 



< 



CO. 



From (8.18) and the above Theorem Proposition 8.1 follows. □ 

Thus we need to check that function (8.11) belongs to N-class, that is, (8.14) 
holds. 

So step 2, will be completed if one proves (8.8), (8.13) and (8.14). 
Assuming (8.8) and (8.13), one proves (8.14) as follows: 

i) First, one checks that (8.14) holds with H(£) in place of h(£). 

ii) Secondly, using (8.8) one writes h(£) as: 



h{£) = H{£) +1 j I [K 1 {r,p)+K 2 {r,p)]u l {p)u t {r)% 



(8.19) 



+ 11 



o Jo 



ds dt 

I 2 ! 2 



Ki(r, t)K 2 (r, s)ue(t)ue(s) . 



Let us now go through i) and ii) in detail. 
In |, 8.411.8] one finds the formula: 



.20) 



7«W^' +1 / (l-t 2 ) e e lrt dt, 



where 7^ is defined in (8.9). 

From (8.20) and (8.10) one gets: 

(8.21) 

»i 



H{t)\ < I dr\p(r)\r 



21+2 



(l-t 2 Ye irt dt 



< ca 2a , £ = a + iT, cr > 0. 



One can assume a > 1 without loss of generality. Note that 
.22) ln + (ab) < ln + a + ln+ b if a, b > 

Thus (8.21) implies 



.23) 



In 



H 



1 - re 1 * 



1 + re} 1 

< |lnc| +21na 

< llncl +21na 



dip < 



In 4 



ca 



d(p 



Re- 



1 - re^ 



1 + re l f 
l-r 2 



dip 



_ 7T 1 + 2r cos ip + r 2 
Here we have used the known formula: 
.24) r dp 2, 



dip < I lnc| + 47rlna < 00. 



< r < 1. 



_ 7r 1 + 2r cos ip + r 2 1 — r 2 
Thus, we have checked that H{£) e N(H + ), that is (8.14) holds for H{£). 
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Consider the first integral, call it in (8.19). One has, using (8.20) and 

(8.13), 

< f dpp-^Kifrp)] + \K 2 (r 7 p)\y+ 1 p t < c{a)a 2 °, 
Jo 

(8.25) £ = a + iT, a > 0. 

Therefore one checks that Ii(£) satisfies (8.14) (with I\ in place of h) as in (8.23). 
The second integral in (8.19), call it I 2 {£), is estimated similarly: one uses (8.20) 
and (8.13) and obtains the following estimate: 

(8.26) \I 2 {£)\ < c{a)a 2a , e = a + ir, a>0. 

Thus, we have proved that h(£) G AT(rj + ). 

To complete the proof one has to derive (8.4) and check (8.13). 

Derivation of (8.4): Subtract from (7.28) with q — qi this equation with q = q 2 
and get: 



1(1+1) ,(2) 
- w - qiw = pip) 



(8.27) W" +W — uj - h1 _uj — iJ y l , 

where 

(8.28) P := qi -q2, w:=i)f\r)-i)f\r). 

Multiply (8.27) by ^^(r), integrate over [0, 00), and then by parts on the left, and 
get 



(8.29) (w'lpP - } ') 



p^^P dr, £ G C. 



By the assumption 5^ = 5^ if I G C, so w and w' vanish at infinity. At r = 
the left-hand side of (8.29) vanishes since 

(8.30) ^ ) (r)=0(r 4+1 ) as r -» 0. 

Thus (8.29) implies (8.4). □ 

Derivation of the representation (8.8) and of the estimate (8.13). One can 
prove |2j| that the kernel K(r,p) of the transformation operator must solve the 
Goursat-type problem 

(8.31) r 2 K rr {r, p) - p 2 K pp {r, p) + [r 2 - r 2 q(r) - p 2 }K{r, p) = 0, < p < r, 

r r 

(8.32) K(r, r) = - sq{s) ds := g(r), 

2 Jo 

(8.33) K(r,0)=0, 

and conversely: the solution to this Goursat-type problem is the kernel of the 
transformation operator (8.8). 

The difficulty in a study of the problem comes from the fact that the coefficients 
in front of the second derivatives degenerate at p = 0, r = 0. 

To overcome this difficulty let us introduce new variables: 

(8.34) £ = ln?' + lnp, r\ = In r — In p. 
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Put 

(8.35) K(r,p):=B(tv)- 
Then (8.31)-(8.33) becomes 

(8.36) Bfr- ^Br, + Q(t,r))B = 0, rj > 0, -oo < £ < oo, 

(8.37) B(£,0) = g(e%) := G(0, 

(8.38) B(-oo,7?) = 0, r?>0, 
where g(r) is defined in (3.32) and 



(8.39) 

Note that 
(8.40) 

(8.41) 



Q(S,v) ■■= 



l 



e s1 ''g e 2 



sup e 2 < c, 

— oo<£<oo 

sup / \Q(s,r))\ds < c(A,B), 

0<V<B J -00 



for any Ael and any B > 0, where c(A, B) > is some constant. 
Let 

(8.42) L(£, rj) ~ B(£,rj)e- * 
Write (8.36)-(8.38) as 

(8.43) L Sn + Q(£, rj)L = 0, 77 > 0, -00 < £ < 00 

(8.44) L(£,0) = e-*G(0:=&(0; £(-00, *?)=<), r? > 0. 

Integrate (8.43) with respect to r\ and use (8.44), and then integrate with respect 
to £ to get: 

(8.45) £ = VL + &, VL := - [ ds f ' dtQ(s,t)L(s,t). 

J -co JO 

Consider (8.45) in the Banch space X of continous function L(£, rj) defined for 
77 > 0, —00 < £ < 00, with the norm 

(8.46) 



\L\\ := \\L\\ab ~ sup (e 7t \L(s,t)\) < 00, 

0<t<B 
— oo<s<A 



where 7 = 7 (A, B) > is chosen so that the operator V is a contraction mapping 
in X. Let us estimate ||V||: 

(8.47) 



\VL\\ < sup 



-oo<£<A J — 00 JO 

0<))<S 



/ rft|Q(s,t)| e - 7(r '-* ) e- 7 *|i(s,t)| 



< ||L|| sup / / f 2e s+t + e s+ * qfe^) ) e^"-*) < -||L||, 

-oo< 5 <aJ-oo Jo v V / / 7 



0<?j<B 



where c > is a constant which depends on A, B, and on f£ r |<?(r) | dr. 
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(8.48) 



If 7 > c then V is a contraction mapping in X and equation (8.45) has a unique 
solution in X for any — oo < A < oo and B > 0. 

Let us now prove that estimate (8.13) holds for the constructed function K(r, p). 
One has 

pT />00 

/ \K(r, p)\p~ 1 dp = r / \L(2 lnr — rj, rj)\e~%dr) < oo 
Jo Jo 

The last inequality follows from the estimate: 

(8.49) \L(Z,V)\ < ce^ +£ ^^+^ i+E2 
where e\ and £2 > are arbitrarily small numbers, 

(8.50) Mi(0:=/ ? <M«)> n(s) := y (l + | 9 (e*) |) 

The proof of Theorem 8.1 is complete when (8.49) is proved. 
Lemma 8.1. Estimate (8.4-9) holds. 
Proof. From (8.45) one gets 

(8.57) m(£,r}) <c + Wm 

where 

(8.52) c := sup \b(0\ < \ [ s(q(s)ds, m(£, rj) := \L(£, v )\, 

-oo<£<oo * JO 

and 

(8.53) Wm:= f 6 ds P dt p(s + t)m(s,t). 

J -co Jo 

Without loss of generality we can take cq = 1 in (8.51): If (8.49) is derived from 

(8.51) with c — 1, it will hold for any c > (with a different c in (8.49)). Thus, 
consider (8.51) with c = 1 and solve this inequality by iterations. 

One has 

ds n(s + t)dt= Hi(£ + t)dt<rip,i(t; + r)). 
-oo Jo Jo 

One can prove by induction that 

(8.55) w^^ mpl . 
Therefore (8.57) with Co = 1 implies 

(8.56) m(^)<l + £ J? "^ + r?) 
Consider 



(8.57) F(z):=l + J2 



00 r 
z 



^ (n!) 2 

n— 1 v ' 



This is an entire function of order i and type 2. 
Thus 



(8.58) \F(z)\ < ce (2+£l)|z| 

From (8.56) and (8.58) estimate (8.49) follows. 



40 



A.G. RAMM 



Lemma 8.1 is proved. □ 
Theorem 8.1 is proved. □ 
9. Discussion of the Newton-Sabatier procedure for recovery of q{r) 

FROM THE FIXED-ENERGY PHASE SHIFTS 

In H and Jl6| the following procedure is proposed for inversion of the fixed- 
energy phase shifts for q(r). We take k = 1 in what follows. 

Step 1 . Given {8e}w=o i 2 one solves an infinite linear algebraic system ((12.2.7) 
in |) 

00 

(9.1) tan<^ = Mw(l + taxi 81 tan 8e')ap 

for constants ag. Here 



if \l — l'\ is even or zero, 

yi 1 r—^ if I* - f I is odd. 



(9.2) Ma> = 

Assuming that (9.1) is solvable and ae are found, one calculates (see formula (12.2.8) 
in |) 

(9.3) 

/ 2 \ I 7ra ^ (l + tan 2 Se) r 1 

c t = ae (1 + tan 8 e ) < 1 - M u ,a e >(ta,ndi> - tan^) > 

I £'=0 J 

Step 2 . If a are found, one solves the equation for K(r, p) (see formula (12.1.12) 

in |) 



(9.4) K(r,s) = f(r,s)- K(r,t)f(t,s)t- 2 dt, 

Jo 

where 

00 

(9.5) f(r, s) := c e u e (r)u e (s), 

e=o 

and ui(r) are defined in (7.22). 

Note that in this section the notations from || are used and by this reason the 
kernel K(r,t) in formulas (9.4) and (9.6) differs by sign from the kernel K(r,p) in 
formula (8.8). This explains the minus sign in formula (9.6). 

Assuming that (9.4) is solvable for all r > 0, one calculates 

/ % , , 2d K(r,r) 

(9.6) q(r) = —- -±L-L. 

r dr r 

Note that if (9.4) is not solvable for some r = r > 0, then the procedure breaks 
down because the potential (9.6) is no longer locally integrable in (0, 00). In |l6) it is 
argued that for sufficiently small a equation (9.4) is uniquely solvable by iterations 
for all < r < a, but no discussion of the global solvability, that is, solvability for 
all r > 0, is given. It is assumed in || and that the sequence {q} in (9.3) does 
not grow fast. In (|j[ (12.2.2)]) the following is assumed: 

00 

(9.7) £|(*|r 2 <oo. 
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Under this assumption, and also under much weaker assumption 
(9.7') \a\ < ce M 

for some b > arbitrary large fixed, one can prove that the kernel (9.5) is an entire 
function of r and s. This follows from the known asymptotics of ui(r) as t — > oo: 

21+1 

r ( er \ 2 1 



(9-8) Mr) = ^- 2 [ WTT ) + 

Thus, equation (9.4) is a Fredholm-type equation with kernel which is an entire 
function of r and s. Since K(r, 0) = and f(r, s) — f(s, r) — at s — 0, equation 
(9.4) is a Fredholm equation in the space of continuous functions C(0,r) for any 
r > 0. 

If (9.4) is uniquely solvable for all r > 0, then one can prove the following: 

Claim. K(r, s) is an analytic function of r and s in a neighborhood A of the 
positive semiaxis (0, oo) on the complex plane of the variables r and s. 

This claim is proved below, at the end of this section. 
Therefore the potential (9.6) has to have the following: 

Property P: q(r) is an analytic function in A with a possible simple pole at r = 0. 

Most of the potentials do not have this property. Therefore, if one takes any 
potential which does not have property P, for example, a compactly supported 
potential q(r), and if it will be possible to carry through the Newton-Sabatier 
procedure, that is, (9.1) will be solvable for ag and generate eg by formula (9.3) 
such that (9.7) or (9.7') hold, and (9.4) will be uniquely solvable for all r > then 
the potential (9.6), which this procedure yields, cannot coincide with the potential 
with which we started. 

An important open question is: assuming that the Newton-Sabatier procedure 
can be carried through, is it true that the reconstructed potential (9.6) generates 
the scattering data, that is, the set of the fixed-energy phase shifts {Si} with which 
we started. 

In ||, pp. 203-205] it is claimed that this is the case. But the arguments in || are 
not convincing. In particular, the author was not able to verify equation (12.3.12) 
in H and in the argument on p. 205 it is not clear why A' t and 5' e satisfy the same 
equation (12.2.5) as the orginal Ag and Si. 

In fact, it is claimed in || (12.5.6)] that S' e = 0(j), while it is known (37| that if 
q(r) = for r > a and q(r) does not change sign in some interval (a — e, a), q =/= 
if r G (a — e, a), where e > is arbitrarily small fixed number, then 

/IP 

(9.9) lim - \8 t \ 

i~*oo v e 

It follows from (9.9) that for the above potentials the phase shifts decay very 
fast as I — > oo, much faster that i. Therefore S' e decaying at the rate j cannot be 
equal to St,, as it is claimed in [|| p. 205], because formula (9.9) implies: 

( ea[l + o (lT 2 

~ { — U 

It is claimed in ||, p. 105, that the Newton-Sabatier procedure leads to "one 
(only one) potential which decreases faster that r~ § " and yields the original phase 
shifts. However, if one starts with a compactly supported integrable potential (or 



42 



A.G. RAMM 



any other rapidly decaying potential which does not have property P and belongs 
to the class of the potentials for which the uniqueness of the solution to inverse 
scattering problem with fixed-energy data is established), then the Newton-Sabatier 
procedure will not lead to this potential as is proved in this section. Therefore, 
either the potential which the Newton-Sabatier procedure yields does not produce 
the original phase shifts, or there are at least two potentials which produce the 
same phase shifts. 

A more detailed analysis of the Newton-Sabatier procedure is given by the author 
in§. 

Proof of the claim. The idea of the proof is to consider r in (9.4) as a parameter 
and to reduce (9.4) to a Fredholm-type equation with constant integration limits 
and kernel depending on the parameter r. Let t = rr, s = ro~, 

, 01fV , K(r, rr) /(rr, ra) /(r, a) 

(9.10) :=b[r;r), := a{a,T,r), := g(a;r) 

T TTO a 

Then (9.4) can be written as: 

(9.11) b(a;r) — g{r;a) — / <z(<7, r; r)6(r; r) dr. 

Jo 

Equation (9.11) is equivalent to (9.4), it is a Fredholm-type equation with kernel 
a(<7, r;r) which is an entire function of a and of the parameter r. The free term 
<?(r, a) is an entire function of r and a. This equation is uniquely solvable for all 
r > by the assumption. Therefore its solution b(o~; r) is an analytic function of r 
in a neighborhood of any point r > 0, and it is an entire function of a [ |35[ . Thus 
K(r,r) — b(l,r) is an analytic function of r in a a neighborhood of the positive 
semiaxis (0, oo). 

□ 



10. Reduction of some inverse problems to an overdetermined 

Cauchy problem 

Consider, for example, the classical problem of finding q(x) from the knowledge 
of two spectra. Let u solve (1.1) on the interval [0,1] and satisfy the boundary 
conditions 

(10.1) u(0) = u(l) = 0, 

and let the corresponding eigenvalues fc^ := A n , n = 1,2, ... , be given. If 

(10.3) u(0)=u'(l)+hu(l) = 0, 

then the corresponding eigenvalues are (j, n , n = 1, 2, . . . . 
The inverse problem (IP10) is: 

Given the two spectra {A ra } U {n n }, n = 1, 2, 3, ... , find q(x) (and h in (10.3)). 
Let us reduce this problem to an overdetermined Cauchy problem. Let 

.s , ,n sin(fca;) f x ,,. s sin (few) , , r ,,. / sinffex) \ 

(10.4) u(x, k) = — i—L + K(x, V)^^- dy ~ (I + K) \— ^ J 

solve (1.1). Then (10.1) and (10.2) imply: 
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and 

(10.6) = cos y/JI^ + K (1,1) 

It is known 0] that 

(10.7) A r , 
and 



Sill y/fl n 



+ / K x (l,y) 



\JThi. 



dy, n = l,2, 



(nn) 2 + c + o(l), n — » oo, 



(10.8) 



7r +ci + o(l), n — > oo, 



where Co and C\ can be calculated explicitly, they are proportional to q(x) dx. 
Therefore, 

(10.9) 



1 + 



tt [ n + 



1 + 



n — > oo. 



It is known jl0| that if (10.9) holds then each of the systems of functions: 
(10.10) {smi^/Kx^n^i^,..., {sin( v / /i^x)} n= i i 2,... 

is complete in L 2 [0, 1]. 

Therefore equation (10.5) determines uniquely {K(l, y)}o< y <i and can be used 
for an efficient numerical procedure for finding K(l,y) given the set {A„} rl= i i 2 j ... 
Note that the system {sin(v / A^"x)} I1= i i 2.... forms a Riesz basis of L 2 [0, 1] since 
the operator I + K, defined by (10.4) is boundedly invertible and the system 
{u(x, \A«)}n=i,2,... forms an orthornormal basis of L 2 (0, 1). 

Equation (10.6) determines uniquely {K x (l,y)}o< y <i if {^ n }n=i,2,... are known. 
Indeed, the argument is the same as above. The constant K(l,l) is uniquely 
determined by the data {fi n }n=i,2,... because, by formula (5.41), 



(10.11) 



1 f 1 

K{1,1) = - J q(x)dx. 



We have arrived at the following overdetermined Cauchy problem: 
Given the Cauchy data 



(10.12) 



{K(l,y),K x (l,y)} Q < y < y 



and the equations (5.40) - (5.4-1), find q{x). 

It is easy to derive [E2L eq. (4.36)] for the unknown vector function 



(10.13) 

the following equation 

(10.14) 

where 



U := 



q(x) 
K(x,y) 



(10.15) 



W(U) 



U = W(U) + h, 

'-2 £ q(s)K(s, 2x - s) ds^ 
| f D q(s)K(s,t) ds dt 
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D xy is the region bounded by the straight lines on the (s, t) plane: s = 1, t — y 
s — x and t — y = — (s — x), and 

7' 



(10.16) h := 

\9 

(10.17) f(x) := 2[K y (l, 2x-l) + K x (l, 2x - 1)], 

K(l,y + x-l)+K(l,y-x + l) 1 , xi 

(10.18) 3 (x):=^^ ' K di L- K s (l,t)dt. 

z z Jy+x-l 

Note that / and g are computable from data (10.12), and (10.14) is a nonlinear 
equation for q(x) and K(x,y). 
Consider the iterative process: 

(10.19) U n+1 = W{U n ) +h, U = h. 
Assume that 

(10.20) q{x) = for x > 1, q = q, ge£°°[0,l]. 
Let xq € (0, 1) and define the space of functions: 

(10.21) L(x ) :— L°°(x Q , 1) x L°°(A Xo ), 
where 

(10.22) A xo :={x,y:xo<x<l, < \y\ < x}. 
Denote 

(10.23) \\u\\ := esssup \q[x)\ + sup \K(x,y)\. 
Let 

(10.24) \\h\\ < R. 

Theorem 10.1 (Q). Let (10.17), (10.18), (10.20), and (10.24) hold, and choose 
any R > R. Then process (10.19) converges in L(xq) at the rate of a geometrical 
progression for any Xq S (1 — /A, 1), where 



fi := mm 

One has 



8(R-R) _2_ 
5.R 2 ' 5.R 



(10.25) lim U n = ( jW),x,yeA X0 . 



K{x,y) / 
If one starts with the data 

(10.26) {K(x Q ,y), K x (x ,y)} < 

replaces in (10.19) h by ho := (g°), where fo and go are calculated by formulas 
(1.17) and (1-18) in which the first argument in K(l,y),x — 1, is replaced by 
x = xq, then the iterative process (10.19) with the new h = ho, in the new space 

L( Xl ) := L°°(xi,xo) x L°°(A Xl ), 

with 

&xx = { x , V : x i < x < Xo, < \y\ < x}, 
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converges to (k(xm)) in L ( x i)- 

Infinite number of steps one can uniquely reconstruct q(x) on [0,1] from the data 
(10.12) using (10.19). 

Proof. First, we prove convergence of the process (10.19) in L(xq). 

The proof makes it clear that this process will converge in L{x\) and that in final 
number of steps one recovers q(x) uniquely on [0,1]. Let B(R) := {U : \\U\\ < R, 

U G L(x )}, R> R. 
Let us start with 

Lemma 10.1. The map U G W(U) + h maps B(R) into itself and is a contraction 
on B(R) ifx G (l-f.,lU:=mm(^,i). 

Proof of the lemma. Let U = ( ) , V = (k 2 ) One has: 

2^(\q 1 -q 2 \\K 1 \ + \q 2 \\K 1 -K 2 \) ds 

hL (l9i-<al|tfil + l<fel l^i-^2|) dsdt 



\\W(U)-W(V)\\ < 
(10.27) <\\U-V\\ 



2(l-x )R+^(l-x ) 2 



< \\U-V\\(l-x )-R. 



Here we have used the estimate (1 — x ) 2 < 1 — x and the assumption \\U\\ < R, 

< R. 
If 



(10.28) 



1 - x < — , 
5R 



then W is a contraction on B(R). 

Let us check that the map T(U) = W{U) + h maps B(R) into itself if 1 — x < \i. 
Using the inequality ab < ^j- if a + b = R, a, b > 0, one gets: 



(10.29) 

Thus if 
(10.30) 

then the map U 
10.1 is proved. 



\\W(U) + h\\ < \\W(U)\\ + \\h\\ 

<2(l-x )^-+ 1 -(l-x f^+R 

^ R 2 
<-(l-x a )—+R<R if 1 



x < 



8(R-R) 
5R 2 



jj, = mm 



8{R-R) 
5R 2 ' 



5R 



W{U) + h is a contraction on B(R) in the space L{xq). Lemma 



□ 



From Lemma 10.1 it follows that process (10.19) converges at the rate of geomet- 
rical progression with common ratio (10.30). The solution to (10.14) is therefore 
unique in L(xq). 

Since for the data h which comes from a potential q G L°°(0, 1) the vector 
sorves (10.14) in L(xq), it follows that this vector satisfies (10.25). Thus, 
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process (10.19) allows one to reconstruct q(x) on the interval from data (10.12), 
xq = 1 — /i, where /i is defined in (10.30). 

If q(x) and K(x, y) are found on the interval (xq, 1), then K(xo, y) and K x (xq, y) 
can be calculated for < \y\ < xq. Now one can repeat the argument for the interval 
(x±, xq), xq — x\ < fj,, and in finite number of the steps recover q(x) on the whole 
interval [0,1]. 

Note that one can use a fixed fi if one chooses R so that (10.24) holds for h defined 
by (10.16) and (10.17) with any x € L°°[0, 1]. Such R does exist if q e L°°[0, 1]. 
Theorem 10.1 is proved. □ 



Remark 10.1. Other inverse problems have been reduced to the overdetermined 
Cauchy problem studied in this section (see |^2| , [p2| . $L% ). The idea of this reduc- 
tion was used in jl2| for a numerical solution of some inverse problems. 

11. Representation of I-function 

The I(k) function (2.1) equals to the Weyl function (2.3). Our aim in this section 
is to derive the following formula ([fL9|): 

J ■ />oo 

(11.1) I{k) =tk + Y, 1^, +"> S: = / a{t)e lkt dt, 

where k :— 0, rj — const > 0, 1 < j < J, r > iff /(0) = 0, a(t) = a(t) is a 
real-valued function, 

(11.2) a{t) S L X (WL + ) if /(0) ^ and q(x) G £i,i(M+), 

(11.3) a(t) S £ X (M+) if /(0) = and q e Li, 3 (K+)- 

We will discuss the inverse problem of finding q(x) given I{k) Vfc > 0. Unique- 
ness of the solution to this problem is proved in Theorem 2.1. Here we discuss a 
reconstruction algorithm and give examples. Formula (11-1) appeared in |is[| . 

Using (2.17) and (2.1) one gets 

(11.4) 

ik - A(0) + L°° Ai(w)e ifey dy 

J W = A r°° ikv a u A (y)'-=A(Q,y), A 1 (y):=A x (0,y). 

1 + Jo e lk yA(y)dy 

The function 



r k 

(11.5) f(k) = l+ A(y)e <k »d V = Mk) 1 -—l[ 

Jo K ■+■ * 



/c i r /c 2 /c j 



/c — 1~ 2 /b-i 

where /o(&) is analytic in C+, /o(°°) = 1 hi C+, that is, 

(11.6) fo(k) -> 1 as |fc| -► oo, and / (fc) 7^ 0, Vfc e C+ := {k : Imk > 0}. 
Let us prove 

Lemma 11.1. If f(0) ^ and q E Li,i(M+) then 

(11.7) f (k) = l+ b Q (ty kt dt:=l + b Q , 6 GlU 14 (M+). 
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Here W /1 ' 1 (R+) is the Sobolev space of functions with the finite norm 

/>oo 

\\b \\w^ := / (\b (t)\ + \b' Q (t)\)dt<oo. 
Jo 

Proof. It is sufficient to prove that, for any 1 < j < J, the function 



)• 



( 1L8 > J ± 7F^ fc ) = 1+ / 9j(t)e ikt dt, 9j eW^ 

K IKj Jq 

Since = 1 + and since A(y) e W 1 ^ 1 (R + ) provided that q e L hl (R + ) 

(see (2.18), (2.19)), it is sufficient to check that 

(11.9) Jt&- = f°° g(t)e lkt dt, geW 1 ' 1 ^). 

K IKj Jq 

One has f{ikj) = 0, thus 



s(k) jm-nik)_ r ... .„.^"'~' t '" - 1 



(11.10) =/ A(y)e- k i y i e^ k - ik ^ s ds= e tks h ] (s)ds 

Jo Jo Jo 

where 

/oo />oo 
Ad/Je-*^"-") = i y A(i + s)e~ k i t dt 

From (ll.ll)one obtains (11.9) since A(y) <= W 1,1 (R+). 

Lemma 11.1 is proved. □ 

Lemma 11.2. 7//(0) = and q E Li j2 (M+), then (11.7) holds. 

Proof. The proof goes as above with one difference : if /(0) = then k = is 
present in formula (11.1) and in formulas (11.10) and (11.11) with ko = one has 

(11.12) ho(s) = i [ A(t + s)dt. 

Jo 

Thus, using (2.18), one gets 



o Jo Jo 



DC 



ho(s)\ds<c / ds dt \q(u)\du 



/•OO />OG />QO />OG />OG 

(11.13) =2c/ (is/ du / \q(u)\du<2c ds \q{u)\udu 

Jo Jf Jv Jo Jf 

/•oo 

= 4c / u 2 |g(u)| cfri < oo if geli. 2 (l + ), 
Jo 

where c > is a constant. Similarly one checks that h' (s) G i 1 (M + ) if q £ 
Li, 2 (K+). 

Lemma 11.2 is proved. □ 
Lemma 11.3. Formula (11.1) holds. 
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Proof. Write 

1 k = l k—ikj 



k+i TT J k+ikj 
A 1 7=1 



(11-14) 777- 

f{k) fo[k) 

Clearly 

(u - i5) k -i 1 u 1 Bw= i +i^ik.' **==<>. *i>°- 

j=l J j=0 J 

By the Wiener-Levy theorem |§|, §17], one has 

( 1L16 ) T77T = 1 +/ Hty kt dt, b(t)£W 1 ' 1 ( 



fo(k) 



Actually, the Wiener-Levy theorem yields 6(i) S i 1 (R+). 

However, since bo £ W i,1 (R+), one can prove that b(t) £ W 1 ' 1 (M.+ ). 

Indeed, b and bo are related by the equation: 

(11.17) (l + 6o)(l + 6) = 1, VfceR, 
which implies 

(11.18) b=~b -b b, 
or 

(11.19) b(t) = -b (t) - f b Q (t - s)b(s) ds := -b Q - b * b, 

Jo 

where * is the convolution operation. 

Since V £ L 1 (R + ) and b £ L 1 (R + ) the convolution b Q * b £ L 1 (M+). So, 
differentiating (11.19) one sees that b' £ L 1 (M. + ), as claimed. 

From (11.16), (11.15) and (11.4) one gets: 

(11.20) 

I(k) = (ik-A(Q)+A 1 )(l + b) I l + y^-Xr- ) =^ + c + 5Zr^V +S < 

I fv iifv'i i rC vivo 

\ 3=0 ] J=0 3 

where c is a constant defined in (11.24) below, the constants a,j are defined in (11.25) 
and the function a is defined in (11.26). We will prove that c = (see (11.28)). 
To derive (11.20), we have used the formula: 

~ikt 



V 6 ^ 

ik 



1 



o ik jo 



DC: 



Akt u' 



e 



b'{t)dt 



= -6(0) - 6', 



(11.21) ikb = ik 
and made the following transformations 

(11.22) I(k) =ik- A(0) - 6(0) -V + M- A(0)b + A^b^2 



k — ik 



^ Cj [A(0) + b(0)} ^ g(k)-Wkj) c g(ikj)c 



3=0 



k %ko k %k n k zkj 

j=0 3 3=0 3 j=0 

where 

(11.23) g(k) := -V + A-y - A(0)b + Aib. 
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Comparing (11.22) and (11.20) one concludes that 

./ 

(11.24) c:=-A(0)-b(0)+i^2cj, 

3=0 

(11.25) cij := -Cj [kj + A(Q) + 6(0) - g{ikj)] , 
(11-26) ~a(k) :=m+t m k ~T i) c - 

3=0 J 

To complete the proof of Lemma 11.3 one has to prove that c = 0, where c is 
defined in (11.24). This is easily seen from the asymptotics of I(k) as k — > oo. 
Namely, one has, as in (11.21): 

(11.27) A(k) = - \a' 

ik ik 

From (11.27) and (11.4) it follows that 

I(k) = (ik-A(0)+A 1 ) l-^+ Q 

(11.28) =(ik-A(0) + A 1 )^l + -^ + o(j^j=ik + o(l), k — > +oo. 

From (11.28) and (11.20) it follows that c = 0. 

Lemma 11.3 is proved. □ 

Lemma 11.4. One has cij = irj, rj > 0, 1 < j < J, and tq = if /(0) ^ 0, and 
ro>0 i//(0)=0. 

Proof. One has 

(11.29) a 3 = Res I(k) 

k=ikj f(ikj) 

From (2.7) and (11.29) one gets: 

(11.30) aj = -^-=i^-:=ir h r, : =^>0, j > 0. 
If j = 0, then 

(11.31) a =Rcs/(fc) : =ft 

fe =° /(0) 

Here by ReSfc = o I{k) we mean the right-hand side of (11.31) since I(k) is, in general, 
not analytic in a disc centered at k = 0, it is analytic in C+ and, in general, cannot 
be continued analytically into C_ . 

Let us assume q(x) £ Li i 2(M+). In this case f(k) is continuously differentiablc 
in C+. 

From the Wronskian formula 

(n.32) mk)f(-k)-.no-k)f(k ) =2i 

k 

taking k — > 0, one gets 

(11.33) /'(0,0)/(0) = -i. 



f'(0,ikj) 
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Therefore if q G L 1 . 1 {R + ) and /(0) = 0, then /(O) 7^ and f(0,0) ^ 0. One 
can prove [jl4[ pp.188-190], that if q G then j^y is bounded as k — > 0, 

From (11.31) and (11.33) it follows that 

* 1 
(11.34) aj = — = tr 0l r := - 



/(0) 



lf(0W 



From (2.17) one gets: 

(11.35) /(0)=*/ A^j/dy. 



Since ^4(j/) is a real-valued function if q(x) is real-valued (this follows from the 
integral equation (5.62), formula (11.35) shows that 

(11.36) [/(0)] 2 <0, 
and (11.34) implies 

(11.37) r > 0. 

Lemma 11.4 is proved. □ 

One may be interested in the properties of function a(t) in (11.1). These can be 
obtained from (11.26), (11.16) and (11.7) as in the proof of Lemmas 11.1 and 11.2. 
In particular (11.2) and (11-3) can be obtained. 

Note that even if q(x) ^ is compactly supported, one cannot claim that a(t) is 
compactly supported. 

This can be proved as follows. 

Assume for simplicity that J = and /(0) ^ 0. Then if a(t) is compactly 
supported then I(k) is an entire function of exponential type. It is proved in [[l7| 
p. 278] that if q(x) ^ is compactly supported, q G L 1 (R+), then f(k) has infinitely 
many zeros in C. The function f'(0,z) ^ if f(z) — 0. Indeed, if f(z) = and 
f'(0,z) — then f(x,z) = by the uniqueness of the solution of the Cauchy 
problem for equation (1.1) with k = z. Since f(x,z) ^ (see (1-3)), one has a 
contradiction, which proves that /'(0, z) ^ if f(z) = 0. Thus I(k) cannot be an 
entire function if q(x) ^ 0, q(x) G L 1 (ffi+) and q(x) is compactly supported. 

Let us consider the following question: 

What are the potentials for which a(t) = in (11.1)? 

In other words, suppose 

,7 

(11.38) I{k) = ik + J2' 



find q(x) corresponding to /-function (11.38), and describe the decay properties of 
q(x) as x — * +00. 

We now show two ways of doing this. 

By definition 

(11.39) f'(0,k)=I(k)f(k), f'(0,-k) = I(-k)f(-k), keU. 

Using (11.39) and (2.23) one gets 

[I(k) - I(-k)]f{k)f{-k) = 2ik, 
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or 

(11.40) 



f(k)f(-k) 



Vfe e 



Iml(k) 

By (2.5), (2.6) and (11.30) one can write (see (l9)) the spectral function corre- 
sponding the /-function (11.38) (y/X = k): 



(11.41) 



dp{\) = { J 



A > 0, 



J2 J=1 2fcjT^(A + fc 2 ) dX, A < 0, 

where 5(A) is the delta-function. 

Knowing dp(\) one can recover q(x) algorithmically by the scheme (5.26). 
Consider an example. Suppose /(0) ^ 0, J = 1, 



(11.42) I{k)=ik + 



k — ik\ 



ik - 



ir\(k + ik\) 
k 2 + k 2 



i[k + 



r\k 



k 2 + k 2 J k 2 + k 2 ' 



Then (11.41) yields: 

(11.43) dp{X) = 

Thus (5.27) yields: 

(11.44) L(x,y) = - 



'f(VX + n^), A>0, 
2k iri S{\ + k{) d\, A < 0. 



rxvA sin VXx sin V~Xy n , shikix) sh(kiy) 
d\- — — — h 2k\T\ - 



A + fc 2 va Va 



and, setting X — k 2 and taking for simplicity 2k\T\ = 1, one finds: 



2n 



dfcfc 2 sin(fc:c) sin(fcy) 



o 



(11.45) 



fc 2 + fc 2 A; 2 
2ri f°° dksm(kx)sin(ky) 

ri f°° cifc [cos k(x — y) — cos fc(ai + y)] 



fc 2 + fc 2 



2fci V e 



-3/1 _ g-fei(a;+l/) 



where the known formula was used 
1 



(11.46) 
Thus 
(11.47) 



cos kx 1 I, 

-5 ^dk = — e am , 

k 2 + a 2 2a 



a > 0, 



fci > 0, 



x G 



L(x,y) 



n 

2ki 



g — kl|as— y| _ e -k t (x+y) 



sh(kix) sh(k\y) 



ki 



Equation (5.30) with kernel (11.47) is not an integral equation with degenerate 
kernel: 

- e -fei|i-3/| _ e -ki(*+») 



(11.48) X(x,y)+ / K(x,t) 



sh{k\t) sh(k\y) 



e -kx\x-y\ _ g-fci (x+y) 



2k\jr\ k\ 
sh(k\x) sh(kiy) 



ki 



dt 



2k\jr\ ki k\ 

This equation can be solved analytically Q , but the solution requires space to 
present. Therefore we do not give the theory developed in p3 but give another 
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approach to a study of the properties of q(x) given I(k) of the form (11.42). This 
approach is based on the theory of the Riemann problem 0. 
Equations (11.40) and (11.42) imply 

(11.49) f{k)f{-k) = l^f, vl:=kl + n. 
The function 

(11.50) f (k) := /(A;) £±4^ ^ in C+. 
Write (11.49) as 

fc - ifci . ..fc + ifei k 2 + k\ 

m kTik x h{ - k) k^iki = WT^- 

Thus 

(11.51) = ^ : =^b)" 

The function fo(—k) 7^ in C_, /o(oo) = 1 in C_, so := j^rrfe] i s analytic in C_. 



Consider (11.51) as a Riemann problem. One has 

, k 2 + k 2 1 [°° „ fc 2 + /t? 

Therefore (see problem (11.51) is uniquely solvable. Its solution is: 

11.53 /o fc = -— A h(k) = t—7 ~, 

k + %v\ k — ik\ 

as one can check. 
Thus, by (11.50), 



(11.54) f(k) = 



k — ik\ 
k + iv\ 



The corresponding S'-matrix is: 
Thus 

1 Z 100 

(11.56) F s (x):=— [1- S(k)}e lkx dk = 0(e- kix ) for as > 0, 

^ J — oo 

F d (x) = Sl e- k ^, 

and 

(11.57) F(x) = F s (x) + F d (x) = O (e"^) . 
Equation (5.50) implies A{x, x) — O (e~ 2feia: ), so 

(11.58) q(x) = O (e- 2kix ) , a ^+00. 

Thus, if /(0) 7^ and a(t) = then q(x) decays exponentially at the rate 

determined by the number k\ , fci = min k; . 

i<j<j 3 
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If /(0) = 0, J = 0, and a(t) = 0, then 
(11.59) I{k) = ik+ l -j-, 



(11-60) f(k).f(-k) = 1 J^—, r > 0. 

Let /o(fc) = Mffl. Then equation (11.60) implies: 
(11-61) Mk)fo(-k) = ^ 2 := r 



and /o(fc) 7^ in C+. 

Thus, since mdg ^j^ = 0, /o(fc) is uniquely determined by the Riemann prob- 
lem (11.61). 



One has: 



/O(*0 - T^, /0(-fc) = T^-, 



and 



(a:) = — / 1 - 1 — e' fex rffc = — / — = 2v e "° x , x > 0, 

and F d (a:) = 0. 
So one gets: 

(11.63) F(x) = F s (x) = 2v e-" ox , x > 0. 
Equation (5.50) yields: 

POO 

(11.64) A(x,y) + 2v A(x, t)e- v ° (t+v) dt = -2v e- v ° (x+v \ y>x>0. 

J X 

Solving (11.64) yields: 

(11.65) A(x,y) = -2"oe-^ x+ » \ + *_ 2voX . 

The corresponding potential (5.51) is 

(11.66) q(x) = O (e- 2 " 0X ) , x^oo. 

If q(x) — 0(e~ fe:E ), k > 0, then a(t) in (11.1) decays exponentially. Indeed, in 
this case b'(t), Ai(y), b(t), A\ * b decay expenentially, so, by (11.23), g(t) decays 
exponentially, and, by (11.26), the function ^^If/j^ '■— h with h{t) decaying 
exponentially. We leave the details to the reader. 
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12. Algorithms for finding q(x) from I(k) 

One algorithm, discussed in section 11, is based on finding the spectral function 
p(A) from I(k) by formula (11.41) and then finding q(x) by the method (5.26). 

The second algorithm is based on finding the scattering data (2.10) and then 
finding q(x) by the method (5.49). 

In both cases one has to find kj, 1 < j < J, and the number J. In the second 
method one has to find f(k) and Sj also, and S(k) = ^feff ■ 

If kj and f(k) are found then Sj can be found from (2.12). Indeed, by (11.1) 

(12.1) % rj := Res I(k) = Z^M. 
From (12.1) and (2.12) one finds 

, ^ ~ „ , 2^ -i 2 o 

(12.2) ' - 



^•[/(ifc,-)] 2 ^-[/(ifc,-)] 2 
If fcj are found, then one can find /(fc) from 7(fc) as follows. Since f'(0,k) — 
f(k)I(k), equation (2.23) implies equation (11.40): 

(12.3) f(k)f(-k) 



Iml(k) ■ 



Define 



(12.4) -w-n^g if m<oo, 

j=i ■? 

and 

(12.5) w(k) := -Z— 17 ^L_^2 if /( )=oo. 

One has 7(0) < oo if /(0) ^ and 7(0) = oo if f(0) = 0. Note that if g e 7i, 2 (R+) 
and /(0) = then .f (0,0) 7^ and /(0) ^ 0. 
Define 

< 1M > "<*> : = {§r 

Then h(k) is analytic in C+, h(k) ^ in C+, and h(oo) = 1 in C+, while h(—k) 
has similar properties in C_. Denote h ^ k ^ '■= h_(k). This function is analytic in 

C_, h-(k) ^ in C_ and h-(oo) = 1 in C— Denote ft(fc) := h+(k). 
Write (12.3) as the Riemann problem: 

(12.7) h+(k) = g(k)h-(k), 
where 

(12.8) ^ = l^k) if AOXoo, 
and 

(12.9) ff(fc ) = 7 _*^^ + l if 7(0) = 00. 

We claim that the function g(fc) is positive for all real k, bounded in a neighborhood 
of k = and has a finite limit at k — even if 7(0) = 0. Indeed, if 7(0) = 0, then 



PROPERTY C FOR ODE AND APPLICATIONS TO INVERSE PROBLEMS. 



55 



/(0) = and one can prove that the function Im jr^ -p- is bounded. Thus the 
validity of the claim is verified. 

The Riemann problem (12.7) can be solved analytically: ln/i+(fc) — lnh-(k) = 
In g(k) and since h+(k) and h-(k) do not vanish in C+ and C_ respectively, In 
and ln/i_(fe) are analytic in C+ and C_ respectively. Therefore 

,12.10) »<*)— '(si£,7^*)' 

(12.11) h(k) = h+{k) if lmk>0, h(k) = h-(k) if lmk<0, 
and 

(12.12) f(k) =w(k)h(k), lmk>0. 

Finally, let us explain how to find kj and J given I(k). 
From (11.1) it follows that 

(12.13) 



1 /-CO 

— J (I(k) - ik)e" Lkt dk = -J2 r ^ ' " ~ > r ' ' ,K 



3=1 



2 



Taking t — > — oo in (12.13) one can find step by step the numbers ro, ki, n, ki, 
r 2 . . . , rj, kj. If 1(0) < oo, then r = 0. 

13. Remarks. 

13.1. Representation of the products of the solution to (1.1). In this sub- 
section we follow Consider equation (1.1) with q = qj, j = 1, 2. The function 
u(x,y) := ipi(x,k)(p2(y,k) where ipj, j = 1,2, satisfy the first two conditions (1.4), 
solves the problem 



(13.1) 



d 2 



d 2 

W 2 ~ q2iy) 



u(x,y) 

"(0, y) = 0, u x (0, y) = <p 2 (y, k) 



u{x,y), 



(13.2) 

(13.3) u(x,0) = 0, u y (x, 0) = <px(x, k). 

Let us write (13.1) as 



(13.4) 



<9 2 



u(x,y) = [qi(x) - q 2 (2/)]u(a;,y) 



\ dx 2 dy 2 ( 

and use the known D'Alembert's formula to solve (13.3)-(13.4): 



(13.5) 



u(x,y) 



[qi(s) — (i)]u(5, t) dsdt ■ 



1 



x+y 



fi(s) ds, 



x-y 



where D xy is the triangle < t < y, x — y + t<s<x + y — t. 

Function (13.5) satifies (13.3) and (13.1). Equation (13.5) is uniquely solvable 
by iterations: 

(13.6) u(x,y) = ^Tu m (x ) y), u (x,y) := - I 



fi(s,x) ds, 



(13.7) 



u m+ i(x,y) = i 



\qi{s) - q 2 (t)]u m (s,t) dsdt. 
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Note that 

(13.8) u m (x,y) = - w m (x,y,s)ip 1 (s)ds. 

^ Jx-y 

If m = this is clear from (13.6). If it is true for some m > 0, then it is true 
to + 1: 

u m +i(x,y) = - dt ds[q-i(s) - q 2 (t)]- / w m (s,ier)^i(cr) da 

1 JO Jx-y+t 1 Js-t 



1 



-y+t 

y r x +v 



2 Jo 

I r'-c+y 



(13.9) =- / dt dcnpi(a)w m (x,y,t,<j) 



x-y 



daip 1 (a)w m+ i(x 1 y,a), 

I x-y 

where w m and w m+ \ are some functions. 

Thus, by induction, one gets (13.8) for all to, and (13.6) implies 

I rx+y 

(13.10) u(x,y) = - w(x,y,s)<pi(s)ds, 

^ Jx-y 

where 

oo 

(13.11) w(x, y, s) := ^ w m (x, y, s). 

To satisfy (13.2) one has to satisfy the equations: 
= / w(0,y, s)cpi(s, k)ds, 

J-y 

(13.12) <P2(y,k) = -[w(0,y,y)ipi(y) - w(0, y, -y)yi{-y)} 

1 f v 

+ 2 J w x (0 7 y,s)tp 1 (s)ds. 

Formula (13.10) yields 

2 r x +v 

(13.13) (p 1 (x,k)(fi 2 (y,k) = - w(x,y,s)(p 1 (s,k)ds. 

J x — y 

If x = y, then 

1 f 2x 

(13.14) (fi(x,k)if2(x,k) = - / w(x,x,s)ipi(s,k)ds. 

2 Jo 

Therefore, if 

/ h(x)(pi(x, k)(f2(x, k) dx = Vfc > 0, 
Jo 

then 

r 2x 



f-a rZx 

= / h(x) / w(x 1 x 1 s)tpi(s,k)dsdx 
Jo Jo 

r-2a r-a 

= I dsLpi(s,k) / dxh(x)w(x, x, s) Vfc > 0. 

Jo ■/# 
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Since the set {<fii(s, &)}v/c>o is complete in L 2 (0, 2a), it follows that = Jf dxh(x)w(x, x, s) 
for all s e [0, 2a]. Differentiate with respect to s and get 

(13.15) U; (i'i' S ) \ k (^) ~ J s dsh ( x ) w s(x, x, s) = 0. 

From Volterra equation (13.15) it follows h(x) — if the kernel w s (x, x, s)w~ 1 (f,f,s) := 

t(x, s) is summable. From the definition (13.11) of w it follows that if \q(x)\ dx < 
oo V6 > 0, then w s (x,y,s) is summable. The function w(x,y,s) has m summable 
derivatives with respect to x,y and s if q(x) has m — 1 summable derivatives. Thus 
one can derive from (13.15) that h(x) = if w (§, |,s) > for all s e [0,2a]. 

If the boundary conditions at x = are different, for example, ip'j(0,k) — 
ho<£j(0, k) = 0, j = 1, 2, then conditions 

(13.16) ?i x — hou\ x _ Q — 0, u y — /iou|^_ =0, /lo = consi > 
replace the first conditions (13.2) and (13.3). One can normalize (fj(x, k) by setting 

(13.17) <^(0,fc) = l. 
Then 

(13.18) i P ' j (0,k) = h o , 

(13.19) u(0,y) = <p 2 {y, k), u(x, 0) = <pi(x, k), 



(13.20) u x (0, k) = h ip 2 {y, k), u y (x,0) = h <fi(x, k), 
and (13.5) is replaced by 

u(x,y) = \ j [qi(s) - q 2 (t)]u(s,t)dsdt 

1 JD xy 

I rx+y i 

(13.21) +- h <pi(s,k)ds+-[<p 1 (x + y,k) + (p 1 (x-y,k)]. 

^ J x — y ^ 

Note that 

1 Id f x+y 1 f x+y 

(13.22) -[ ( p 1 (x + y) + < pi ( x -y)] = -—j <p 1 (s)ds = -J Vs (s)ds. 

Equation (13.21) is uniquely solvable by iterations, as above, and its solution is 
given by the first formula (13.6) with 

h f x+y 1 f x+y 

(13.23) u (x,y) = — ipi(s,k)ds+ - / ip u (s,k)ds. 

^ J x — y ^ J x — y 

The rest of the argument is as above: one proves existence and uniqueness of 
the solution to equation (13.21) and the analog of formula (13.10): 

r x+y 



2 rx-t-v 

(13.24) u(x,y) = - w{x,y,s)$(s)ds, <& := h a ipi{s, k) + ipu(s, k). 

^ Jx-y 

Thus 

p2x 

(13.25) u(x,x) = <pi(x, k)(p 2 (x, k) = / t(x, s)[ho<fi(s, k) + <pu(s, k)] ds, 

Jo 
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where 

(13.26) t(x, s) :— —w(x, x, s), 

and t(x, s) is summable. 

Thus, as before, completeness of the set of products {ipi(x, k)ip2(x, &)}vfc>o can 
be studied. 

13.2. Characterization of Weyl's solutions. The standard definition of Weyl's 
solution to (1.1) is given by (2.2). 
In |n| it is proved that 

(13.27) W(x,k) = e ifcE (l + o(l)) as \k\ -»■ oo, \x\ < b, k 2 £ A, 
where A := {A : \ImX\ > e, dist(X,S) > e}, 

/>±oo 

(13.28) S':=MU[«7_,«7 + ] ) 7± := inf ±/ [u' 2 + g|u| 2 ] dx. 

u(0)=0 

The relation (13.27) gives a definition of the Weyl solution by its behavior on 
compact sets in the x-space as \k\ — > oo, as opposite to (2.2), where k is fixed, and 
x ~ * oo. For multidimensional Schrodinger equation similar definition was proposed 
in jl7|, p. 356, problem 8]. 

We want to derive (13.27) for potentials in Li i i(R_|_) and for k > 0, k — » +oo. 

The idea is simple. For any q = q 6 L i 1 oc (R+), one can construct ip(x,k) and 
?/>(£, k), the solutions to (1.1) and (1.4), for any |ac| < b, where b > is an arbitrary 
large fixed number, by solving the Volterra equations 

, ,n sin(fca;) f x sin[fc(a; — y)] , . , , . , 

(13.29) ip{x,k) = — Y^ + J q \ y)i q(yMy,k)dy 

(13.30) *) = cos(fcx) + ^ ""^"^ gfoMy, fc) dy. 
One can also write an equation for the Weyl solution W: 

(13.31) W(x, k) = cos(kx) + m(fc)^M + f Sm ^ X ~ y ^ q(y)W(y, k) dy. 

k Jo 

This equation is uniquely solvable by iterations for \x\ < b. 
It is known that 

(13.32) m(k) = ik + o(l), \k\ ->• oo, Jmfc > e|JJefc|, e > 0. 

For g 6 Li,i(R + ) the above formula holds when k > 0, k — ► +oo. From (13.31) and 

(13.32) one gets, assuming fc > 0, 

(13.33) fc) = e fc (l + O (£) ) + £ ^Ife^i g (y)W(y, k) dy. 

Solving (13.33) by iterations yields (13.27) for k > 0, k -> +oo. For q e Li ; i(M + ) 
the Weyl solution is the Jost solution. Therefore the above result for k > 0, k — > +oo 
is just the standard asymptotics for the Jost solution. It would be of interest to 
generalize the above approach to the case of complex k in the region (13.27). □ 
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One can look for an asymptotic representation of the solution to (1.1) for large 
|fc|, Imk > s\Rek\, e > 0, of the following form: 

(13.34) u(x, k) = e tkx+ !« ff(t < fe) dt , 
where 

(13.35) a' + 2ika + a 1 - q(x) = 0, cr = ^ + o ' ' 



2ik \k 

From (13.34) one finds, assuming q(x) continuous at x = 0, 

If q(x) has n derivatives, more terms of the asymptotics can be written (sec [ fl4] , 
p.55]). 

13.3. Representation of the Weyl function via the Green function. The 

Green function of the Dirichlet operator L q — — 4^ +q(x) in L 2 (R + ) can be written 
as: 

(13.37) G(x,y,z)=tp(y,y/Z)W(x,Vz), *>V 

where ip(x, k), k := s/z, solves (1.1) and satisfies the first two conditions (1.4), and 
W(x, y/z) is the Weyl solution (2.2), which satisfies the conditions: 

(13.38) W(0,y/z) = l, W'{Q,«Jz) = m{«Jz). 
From (1.4), (13.37) and (13.38) it follows that: 

ncqn\ d 2 G(x,y,k) 

(13.39) — — = m{k). 

OXOy x=y=0 

If q G Li,i(M + ) then W(x,k) = ^jjj^-, where f{x,k) is the Jost solution (1.3), 
k e C+. Note that (13.17) and (13.38) imply: 



f°° 6{x y t) 

G(x,y,z) = ' ' dp(t), 6(x,y,t) = ip(x,Vi)(p(y,Vi). 



(13.40) -Im G(x, y,X + iO) = — ' >YKy ' ' Im m(v / A + iO) 

7T 7T 

and 
(13.41) 

Thus 

(13.42) -Im G(x, y, t + i0)dt = 9(x, y, t)dp(t). 

7T 

From (13.42) and (13.40) one gets, assuming p(— oo) = 0, 

i ft 

(13.43) p (t) = - Imm(VxTi0)dX. 

n J-oo 

If A < then Imm (VA + iO) = except at the points A = — fc 2 at which 



fiikj) = 0, so that m(\/ — k 2 + iO) = oo. Thus, if t and a are continuity points of 
p(t), then 



1 '* 



(13.44) p(t)-p(a) = - Imm{V\ + i0) dX, a > 0. 

I" J a 

Let us recall the Stieltjes inversion formula: 



60 



A.G. RAMM 



If z — a + ir, t > 0, p(t) is a function of bounded variation on '. 

(13.45, *)=-/"f*a. 

and if a and b are continuity points of p(t), then 

(13.46) -/ Imtp(X + iO)d\ = p(b) - p(a). 
Therefore (13.44) implies 

(13.47) m( V5) - r ^ 



oo * ' 



The spectral function dp(t) does not have a bounded variation globally, on the 
whole real axis, and integral (13.47) diverges in the classical sense. We want to 
reduce it to a convergent integral by subtracting the classically divergent part of it. 

If q(x) = 0, then p := p (t) for t < 0, m(VX) = iVX, and formula (13.44) with 
a = yields 

(13.48) Po (X) = 

If q(x) = then G(x, y, X) = sin <£ y) e^ x , y < x, so (13.39) yields m{s/X) = 
i>/X. Formula (13.47) yields formally 

(13.49) iVX = - / - — -. 

7T J Q t — X 

This integral diverges from the classical point of view. Let us interpret (13.49) as 
follows. Let ImX > 0. Differentiate (13.49) formally and get 

(13.50) _J_ = If°_^_ ImX>0. 
2s/X Wo (i-A) 2 ' 

This is an identity, so (13.49) can be interpreted as an integral from to A of 
(13.50). The integral J °° t~ * dt which one obtains in the process of integration, 
is interpreted as zero, as an integral of a hyperfunction or Hadamard finite part 
integral. 

Subtract from (13.47) the divergent part (13.49) and get: 
(13.47') m{Vz)-iVZ= r 

J -co t - z 

where 

(13.47") rfcr(A) = dp(X) - dpo(X), dp (X) := 

Integral (13.47') converges in the classical sense if q e Li ) i(R+). Indeed, by (2.5) 
and (13.47") one has da(t) = ^. - l)dt. By (5.65) one has f{Vt) = 1 + 

O(^) as t -> +oo. Thus efcr(i) = 0(-^)dt as t -» +oo. Therefore integral (13.47') 
converges in the classical sense, absolutely, if Imz ^ 0, otherwise it converges in 
the sense of the Cauchy principal value. 
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Let us write (11.1) as 
(13.51) 



m(k) — ik 



iki 



■J2 rj e^H(-t) + a(t)H(t) 
From (13.51) and (13.47') one gets 



(13.52) 
where 
(13.53) 



r Ms) = r 

J-oo S-\ J-, 



e lkt a{t) dt, 



dt, H(t) 



■iO, 



1, t>0, 
0, t<0. 



«(*) :=-S^ J ^(-t) + «(W). 
j=o 

Taking the inverse Fourier transform of (13.52) one can find a(t) in terms of 
a(s). If k > then k = y/k 2 + iO and if k < then fc = Vfc 2 - iO. Thus: 

cfcr(s) 



1 z 100 r 



s — k 2 — iO 



(13.54) 



1 

"27T 



>'/° 

Jo 



ifci 



dfc- 



+ 



/" 

J — oc 



dk- 



k 2 + iO — s ./_™ k 2 — iQ — ,s- ' 

Let us calculate the interior integral in the right-hand side of the above formula. 
One has to consider two cases: s > and s < 0. Assume first that s > 0. Then 
e -ifct /-o 



/■oo g— ifct /"0 p — ifct /• 

J k 2 + i0-s J^ k 2 -i0-s J_ 



-ikt 



dk 



k 2 - 



I 

Jo 



e~ lkt 8{k 2 -s)dk- 



(13.55) 
where 



I 



e- tkt S(k 2 - s)dk 



J :-- 



— F sm.(y/st)H{s) + J, 
V s 

~ ikt dk 



k 2 



1^1 1*1 



J = F sin(|i|-v/i). 

V s 



If s < 0, then 
(13.56) 

If s > 0, then 
(13.57) 

From (13.54)-(13.57) one gets 

/■OO 

(13.58) a(t) = / dcr(s) 

Jo 

Formula (13.58) agrees with (13.53): the second integral in (13.58) for t > is 
an i 1 (R + ) function, while for t < it reduces to the sum in (13.53) because 
dcr(s) = dp(s) for s < 0, dp(s) for s < is given by formula (2.5) and the relation 
between Cj and rj is given by formula (2.7). 



sin(yi) gw iy° do(s) fi _ |t| ^^ 
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